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(KN)  theory  to  predict  absorption  and  scattering  coefficients  of  powders  from  val-j 
lues  of  diffuse  reflectance  of  an  Incident  radiation  beam  is  shown  to  lead  to  er- 
rors of  20%  or  less,  for  realistic  model  powders.  The  use  of  a conventional  plaus-] 
ible  relation  to  predict  imaginary  refractive  index  from  absorption  coefficient  is 
found  to  produce  errors  of  the  order  of  200%,  Mctlnds  for  correcting  these  errors 
are  discussed.  A new  radiative  transfer  equation  is  derived  and  analyzed  for  ap- 
(plicatlons  to  close-packed  strongly  absorbing  or  strongly  reflecting  media.  For 
such  media,  the  new  equation  predicts  reflectances  and  transmittanc.es  which  are  - 


DD 


FORM 
AN  7] 


1473  rOlTlON  OF  1 NOV  65  IS  OBSOLETE 

''^0X00. 


Unr*  i riea 

SECURITYCL  ASsTf  I C 1 O N of  this  PAUt  D*r«  tlntfmt)  ^ 


tCCUWiTY  CLASSIFICATION  OF  TH:S  PAOE(TWi.o  £)•<•  £nlmr»<S) 


Block  20.  (continued) 

significantly  different  from  those  predicted  by  the  standard  transfer  equation. 
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I.  Introduction 

The  major  objective  of  the  theoretical  investigation  reported  here 
was  to  obtain  relations  between  the  actual  absorption  and  scattering, 
coefficients  of  powders,  and  tlie  predictions  of  tlicse  coefficients  which 
are  often  made  by  applying  the  Kubelka-ilunk  (K>I)  theory^  to  measured 
values  of  diffuse  reflectance  from  powder  samples.  Another  objective  was 
> evaluate  >)ne  of  the  conventional  nietho.is  far  obtaining  i.naginarv  refrac- 
tive indices  from  absorption  coefficients,  and  to  suggest  and  implement 
■•orreccions  to  this  metaod,  if  possible.  A third  objective  was  to  asc 'r- 
tain  to  what  extent  the  standard  Mie  theory  and  radiative  transfer  theory 
predictions  of  diffuse  reflectances  from  model  powder  samples  agree  with 
experimental  results  from  corresponding  actual  powder  samples. 

The  first  and  second  objectives  above  were  partially  attained.  The 
results  of  the  investigation  show  that  the  KM  theory  predicts  values  of 
absorption  and  scattering  coefficients  that  are  sensitive  to  incident 
beam  direction  and  sample  optical  depth,  and  that  the  abovementioned 
method  for  conversion  from  absorption  coefficients  to  imaginary  refractive 
Indices  is  subject  to  large  errors.  However,  the  third  objective  mentioned 
above  was  not  achieved,  inasmuch  as  actual  experimental  results  from  powders 
which  could  be  modelled  easily  were  not  available. 

Tlie  investigation  proceeded  on  two  levels.  The  first  level  involved 
computer  simulation  of  several  realistic  powder  samples,  using  standard 
Mie  theory  and  radiative  transfer  theory  algorithms,  wiilch  allowed 
numerical  comparisons  of  predicted  vs.  actual  absorption  and  scattering 
coefficients,  and  evaluatfon  of  the  conventional  relation  between  imaginary’ 

I 

index  and  absorption  coefficient.  The  second  level  involved  the  derivation 
and  fundamental  study  of  a new  equation  of  radiative  transfer,  which  should 
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be  significantly  more  accurate  than  the  standard  one  in  practical  applica- 
tion to  closely  packed  random  scattering  media  which  contain  strongly 
hsorhliig  or  strongly  reflecting  particles. 

lliis  report  is  structured  for  readers  of  three  different  levels  of 
interest.  The  reader  who  needs  only  a synopsis  just  a bit  more  detailed 
the  abstract  n-:-d  consider  only  section  IT,  with  perhios  a reference 
to  the  glossary  of  symbols  in  section  III.  ITie  reader  t>;ho  wants  a look 
at  . rie  ccuiput  ational  'results,  or  it  the  theur-ticii  aiiiceic.s  leading  to  the 
new  equation  of  radiative  transfer,  can  peruse  sections  IV  and  V.  Finally, 
the  reader  who  is  interested  in  depth  can  delve  into  the  appendixes. 

Some  of  the  material  in  this  report  will  soon  be  submitted  for 
publication  in  J.  Opt.  Soc.  Am.  and/or  Applied  Optics.  It  is  planned  to 
write  two  papers  based  on  this  investigation,  one  dealing  with  the  new 
radiative  transfer  equations,  the  other  with  the  computational  results. 

II.  Summary  of  Results  and  Prognosis 
A.  Summary 

The  major  results  of  this  theoretical  investigation  are  as  follows: 

1)  The  use  of  the  Kubelka-Munk  (KM)  theory  to  obtain  values  of  the 

ab.sorption  and  back-scattering  coefficients  of  random  scattering  media 

from  diffuse  reflectance  data  leads  to  errors  < 20%,  over  a range  of 

incident  beam  angles.  These  errors  result  primarily  because  the  KM 

theory  neglcctr.  the  effect  of  incident  beam  angle.  Tlie  errors  might  be 

2 

partially  removed  by  making  use  of  the  Reichman  theory.  See  section  IV 
and  appendix  C for  details. 

2)  The  use  of  the  plausible  hypothetical  relation  n = (Nv) (47rK/A) , 
in  order  to  determine  a value  of  the  imaginary  index  k of  a substance 
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from  a presumed  known  value  of  the  absorption  coefficient  a of  a random 
scattering  medium  composed  of  particles  of  that  substance,  is  subject  to 
errors  of  a factor  of  two  to  tlirec.  In  most  cases,  these  large  errors  can 
be  removed  only  by  recourse  to  the  full  Mie  theory  relation  between  a 
and  <,  which  is  very  difficult  to  invert  in  general.  See  section  TV  and 
appendix  D for  details. 

3)  The  use  of  the  standard  differential  equation  of  radiative  transfer 
for  closely  packed  strongly  reflecting  or  absorbing  scattering  media 
probably  leads  to  significant  error.  For  such  media,  the  new  difference 
equation  of  radiative  transfer  which  was  developed  during  this  investiga- 
tion should  be  used.  The  use  of  the  in  principle  incorrect  far  field 
phase  function,  rather  than  a more  nearly  correct  near  field  one, 
probably  makes  no  pratical  difference.  See  section  V and  appendixes  E, 

A,  1 for  details. 

B.  Prognosis 

Future  work  in  this  area  should  Include  the  following: 

1)  A thorough  Investigation  of  the  new  difference  equation  of 
radiative  transfer,  and  comparison  with  the  standard  theory,  for  realistic 
c]o<-H-parh*‘d  model  seal  ! eri  ng  medj.a  of  prartic.il  inKrest.  Sner  medir. 
include  both  highly  reflecting  and  highly  absorbing  paints  and  coatings, 
powcers  rontaining  sfonp  reflectors  or  ahserbers,  anc  th.  lie.  . 

2)  Another  look  at  the  conventional  averaging  over  different  kinds 
and  sizes  of  particles  in  close-packed  media  composed  of  strongly  absorb- 
ing or  strongly  reflecting  particles.  Prcliminarv  consideration  bv  the 
author  indicates  that  the  standard  method  of  averaging  (see  appendix  A) 
may  lead  to  significant  errors,  for  such  media,  in  addition,  the  effects 
of  polarized  light  and  non-spherlcal  scatterers  should  bo  evaluated. 
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3)  An  error  analysis  for  the  application  of  tlie  Reichman  theory 
(appendix  C)  , rather  than  the  KM  theory,  in  reduction  of  -experiTiiental 
diffuse  reflectance  data,  in  cases  where  the  incident  beam  an^le  is 
not  near  60“ . 

4)  An  attempt  to  invert  tlie  Mie  tlieory  calculations,  so  that  an 
imaginary  and/or  real  index  of  refraction  could  he  obtained  easily  and 
accurately  from  known  values  of  the  scattering  and  absorhtion  cross- 
sections.  (See  appendix  D) . 

III.  Glossary  of  Symbols 
X ^ wavelength,  pm 

(ni,<^)5  (real,  imaginary)  refractive  indices  of  i*^^  optical  type  in  a 
random  scattering  medium 

m^  “ n^  - i<^  = complex  refractive  index  of  particle  type 
fj(r)  = size  distribution  of  spherical  particles  of  type  1 as  function 

N H total  number  density  in  medium  = 

3 ^ til  3 

v^  H dr  f^(r)(4tir  /3)=  average  volume  of  1 type  particle,  nm 

o(f)  c!  i f f f•r^■n  t ) a]  scettcrinj,  ' ross-sec  ! ion  cl  ti-.c  c;vi  •■ay',.'  pa  ' ; j i , i 

in  a meditim  as  function  of  scattering  angle  0 

(s,s')^  di  f i f renti.-.J  crc'ss-sect  i on  lor  scattering  frov,  rlirecticn  ’ to 

direction  s. 


of  particle  radius  r;  dr  f.(r)  = 1 

J th  -3 

number  density  of  particles  of  i type,  pm 


^^abs’°sca’°ext^  ^(absorption,  scattering,  extinction)  cross-sections  of 


2 

an  avf.r;:g(  particle  in  the  medium,  (imi) 

(a,P,Y)  ^(absorption,  scattering,  extinction)  coefficient  of  the  medium. 
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(o) 


(j^  B/y  -■  sinj’] e-scatrering  albedo  of  the  medium. 
p(G)  = phase  function  of  the  medium.  p(9)  3 47ra(0)/a^^^ 

(p,p')  azimuth  averajjed  phase  function  of  the  mediiim,  as  function 
of  )J  = cosO,  p'  = cosO',  for  scattering  from  direction  0' 
to  direction  0.  Here,  0,  O'  are  measured  from  the  forward- 

direction  normal  to  the  plane  surface  of  a sample. 

^ 1 I (o)  / , IN  .forward. 

dp  p (p,±p  ) = average  ( ^ ) -scattering 

o ^ 


- 2 ^ 


dp 


coefficient,  pm 

3 

D ::  mass  density  of  model  powder,  gm/cm  . 

3 

p - mass  density  of  type  i powder  constituent,  gra/cm 

0 H angle  of  incident  beam  w'.r.  to  the  forward-direction  plane 
surface  normal  to  the  sample. 

(K,S)  = Kubelka-Munk  (KM)  (absorption,  scattering)  coefficients, 
pm~^.  K = 2a,  S H 2g-. 

d = Tliickness  of  sample;  optical  thickness  t = yd* 

(K,S)  r KM  (absorption,  scattering)  coefficient  evaluated  from  diffuse 
reflectance  data,  for  each  (0^,d)  considered. 

= diffuse  reflectance  of  a plane  parallel  sample  thickne'^s  d, 
on  ;i  '.'•ark  backcround  . 

R “ diffuse  rc'flectance  of  an  infinitilv  thick  sample. 

TV.  Computer  Modelling  of  Powders 
A.  Method 

Given  the  number  densities,  optical  constants,  and  size  distributions 

ot  ilie  splierical  particles  in  a random  scattering  medium,  standard  Mie 
3 

theory  algoritiims  yield  tlie  quantities  o , , a , o , a,B,Y.o(G)  and 

abs  sea  ext  ’ » i ' 

p(0).  However,  for  closely  packed  media,  such  as  powders,  a given  particle 
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I 

sees  the  near  zone  or  induction  zone  fields,  rather  than  the  radiation 
1 i zone  fields,  from  its  neighbors.  Therefore,  it  is  likely  that  a suitably 

. defined  near  field  phase  function  should  ho  used;  see  appendix  A for 

details.  In  this  work  calculations  were  done  with  both  the  standard 
I asymptotic  phase  function,  and  a near  field  phase  function,  and  compared. 

For  an  actual  powder  sample,  the  particle  number  densities  N.  are  not 
* known  a^  pr i (Vri , but  the  total  mass  density  and  the  fraction  by  weight 

I A of  each  constituent  can  be  regarded  as  given.  Appendix  F details  the 

, conversion  of  these  data  to  N. . 

I I In  an  actual  powder,  the  size  distributions  of  the  various  components 

are  difficult  to  determine  accurately.  Only  good  estimates  of  the  mean 

I 

I particle  radius  and  the  width  of  each  distribution  may  be  available.  In 

modelling  powders,  therefore,  recourse  is  made  to  several  size  distributions 

;j  which  are  thought  to  be  realistic.  'The  various  distributions  used  in  this 

fc  work  are  detailed  in  appendix  G;  these  distributions  have  been  used  recently 

I 4 

j by  other  authors.  It  is  important  to  remark  here  that,  for  the  computer 

I I 

I q modelling  comparisons  which  are  made  in  this  work,  the  actual  size  distribu- 

' _ tions  are  irrelevant,  since  no  comparisons  are  made  with  actual  experimental 

fiat.';.  bf'wevf  r.  rea  1 i • t i r ent  ir<  s were  matir  1 'ir  flu  d:vt’'ii'rt  ii'n  nararu  t • rt 
A in  the  modelling. 

’•'roTv.  tiiij  phase  iar.ction  p(‘- ),  tin  azimvilh  averaged  phase  f\;nction 

' I ^(h,h')  is  found  by  the  method  detailed  in  appendix  H.  Since  only  total 

* / \ 

■ forward  and  backward  radiation  fluxes  are  needed  in  this  work,  only  p °'^(y,u') 

■ is  needed,  i at  li<  r t li.in  t hi>  J itl  1 .tzimiuh-depcndt  nt  ]'iia,'.:(;  frnctien. 

C Given  p (p,p'),  the  thickness  d of  a plane  parallel  sample,  and  the 

Incident  radiation  beam  direction  G , standard  radiative  transfer  theory 
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algorithms  yield  the  total  diffuse  reflectance  and  transmittance  of  the 
sample.  T1)C  doubling  method  algoritlims  which  were  used  in  this  work  are 
detailed  in  appendix  J;  tlie  relevant  radiative  transfer  equation  to  which 
these  were  applied  is  given  in  appendix  A,  Eq . (A19).  A method  used  by 
experimentalists^’^  to  obtain  values  of  the  absorption  and  scattering 
coefficients  of  a medium  from  diffuse  reflectance  data  involves  the  Kubelka- 
Munk  (KM)  theory.^  Tliis  theory  is  detailed  in  appendix  B.  In  this  investi- 
gation, the  diffuse  reflectances  (R^.R^),  wliich  were  calculated  as  described 
above,  were  inserted  into  the  KM  relations  (B2) . The  resulting  quantities 
K(0^,d),  S(0^,d)  were  compared  with  the  knovm  (calculated)  values  of 
K 2 2a,  S = 20_,  for  each  model  powder  co..sidcred. 

A two-flux  radiative  transfer  theory  which  takes  incident  beam  direc- 

2 

tion  into  account  has  been  developed  by  Reichman.  This  theory  is  described 

in  appendix  C.  In  this  invetigation,  the  accuracy  of  this  theory  was  checked 

by  comparing  its  predictions  of  R^  withthoseof  the  accurate  doubling  method. 

As  discussed  in  appendix  C,  it  appears  that  experimentalists  might  achieve 

better  approximations  to  (K,S)  by  using  measured  values  of  (R,,R  ) with 

d oo 

this  Keichman  theory  than  they  can  with  the  standard  KM  the-rv.  i ; thov 

arc  limMcU  t (;  an  incident  heat"  direction  cuitf  riific'rnt  i rt'r  fcl,”. 

o 

The  computer  programr-  v]  Irh  v’ere  developed  for  this  vorV  art  dt  scribed 
hrielly  in  appendix  K.  Complete  documentation  and  listing  is  available 
from  the  author  on  request. 

B.  Re.nults 

Numerical  results  arc  reported  bore  ft>r  a selected  number  of  samples 


and  wavelengths. 
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The  three  saraples  reported  were  constituted  as  follows:  Sarapl e 1, 

BaSO^;  S.unple  2,  BaSO^  plus  carbon;  Sample  3,  model  dust  + BaSO^;  that 

is,  the  model  dust  contained  carbon,  anaionium  sulfate,  and  clay  minerals. 

Table  1 gives  the  relevant  parameter  values  wliich  were  assigned.  lliese 
8 

values  were  quoted  as  typical  for  these  material  constituents  of 
atmospheric  dust  aerosols.  The  BaSO,  was  included  in  all  these  model 
samples  so  tliat  they  would  be  realistic  models  of  actual  powder  samples 
used  in  some  diffuse  reflectance  spectroscopy  experiments. 

The  BaSO^  particles  have  a fairly  narrow  size  distribution  centered 
around  0.5  pm  mean  radius,  approximately.  The  carbon  and  (NH^)2S0^  are 
similarly  sharply  peaked,  but  arc  centered  roughly  around  0.25  pm 
radius.  The  clay  minerals  distribution  is  fairly  broad,  between  0.03  pm 
and  5.0  pm. 

The  values  used  for  the  refractive  indices  n,  ic , and  the  mass  densities 
p,  are  in  agreement  with  those  given  in  several  handbooks,  such  as  that  of 

9 

the  Geological  Society  of  America  (GSA) . However,  the  GSA  handbook  lists 

3 

the  mass  density  of  graphite  as  2.267  gra/cm  , significantly  larger  than 

3 

the  figure  1.7  gm/cm  which  was  used  for  c.arbon  soot  in  this  work.  Tt  v’.as 
III  id((-  ;i  (■>  the  1 t'wr ' value,  l;h.  one  quol.g  to  ;tr.;hor  ;i.'-  art  ropr  i a ; 1 1 

for  soot.  Ag.iin.  it  is  worth  emphasizing  that  the  t rercks  revealed  bv  the 
following  results  arc  insensitive  to  the  accurticy  of  t tu  mooellir.';  of  an 
actual  powder  sample. 

The  following  tables  display  the  results  of  the  computer  calculations 
for  1 1)C5  e model  samples,  at  w.^vel  engths  A = O.B  pm  .in'-  1.0  ,m. 
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Table  A.  Total  diffuse  reflectances  for  the  mocT^l  samples  from 
tlie  accurate  multiple  scattorinp,  calcu'  itions  • ieichman  theory 

(appendix  C) . 

Wavelength  A = 0.6  pm.  Tar  field  phase  functions  were  used. 


Sample 

Components 

T 

0 

Kof  Icctances  , 
cl 

Multiple 

Scattc^ring 

Reichman 

1 

BaSO^^ 

0“ 

0.7306 

0.82)7 

16 

12° 

0.7309 

0.8234 

60° 

0.8155 

0.8595 

0° 

0.9184 

0.9503 

64 

12° 

0.9185 

0.9507 

o 

0 

0.9441 

0.9608 

0° 

0.9773 

0.9871 

256 

12° 

0.9774 

0.9873 

60° 

0.9845 

0.9899 

2 

BaSn^^ 

0° 

0.7101 

0.7976 

Carbon 

16 

12° 

0.7104 

0.7995 

60° 

0.7982 

0.8397 

0° 

0.8406 

0.8745 

64 

12° 

0.8407 

0.8756 

60” 

0.8872 

0.9003 

0° 

0.8460 

0.8/58 

256 

12° 

o 

CO 

ro 

0.8769 

60° 

0 . 8909 

0.9013 

3 

i5aS0, 

0° 

0.7280 

0.8186 

Dust 

16 

1 2° 

0.7283 

0.8203 

60° 

0.8132 

0.8570 

0° 

0.9057 

0.9370 

64 

12° 

0.9058 

0.9375 

60" 

0.9349 

0.9502 

0° 

0.9400 

0.9533 

256 

12° 

0.9401 

0.9537 

60° 

0.9584 

0.9631 

I 


e 

;i 

CM 

• 

O 

c 

f-H 

o 

‘O 

X) 

• 

u 

o 

cO 

O 

II 

(0 

O 

i-H 

r. 

-d-6 

o 

;t 

C/J 

t;3 

o 

« 

in 

CO 

t-l 

CM 

m 

<U 

o 

iM 

CL 

II 

K 

oj 

cn 

CO 

U‘ 

o 

1 

E 

fH|c>J  r< 

m M x: 

u 

<y  -ti  to 
t-^  <u  c 

Xi  C) 

(TJ  •>  rH 

H <u  a 

u > 

0)  eg 

PC  pi 


I 

I 


Tablr  6.  Noar  fir-ld  vs.  far  fiald  comparisons  of  reflectances 
p.'ven  by  the  accur.it.e  inultlple  scattering  calculation  and  by  the 
I’  ■ chman  tlicory  for  siimplc  2 (TiaSO^  plus  Carbon). 

Wavelength  X = 0.6  pm. 
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Tabic  9.  Total  diffuse  reflectances  for  the  nK>del  samples, 
from  the  accurate  multiple  scattering  calculations  vs.  the  Keichinan 
Theory  (appendix  C) . 

Wavelength  A = 1.0  pin. 


Sample  Op! 

T £ 

;cal  depth 
yd 

0 

o 

Keflectances 
Multiple  Rcicfiman 

Scatter ing 

1 

0" 

1 

0.75/0 

0.8368 

BaSO. 

16 

3 2“ 

0.7603 

0.8384 

60“ 

0.835i 

0.8734 

0“ 

0.9280 

0.9550 

64 

12° 

0.9290 

0.9555 

60“ 

0.9512 

0.9651 

0“ 

0.9800 

0.9884 

256 

12“ 

0.9803 

0.9885 

60“ 

0.0864 

0.9910 

2 

0“ 

0.7381 

0.8145 

BaSO, 

16 

32“ 

0.7415 

0.8163 

and 

60“ 

0.8198 

0.8555 

Carbon 

0“ 

0.8572 

0,8855 

64 

12“ 

0.8590 

0.8866 

60“ 

0.9002 

0.91 06 

0“ 

0.86J9 

0.8868 

256 

12“ 

0.8636 

0.8878 

60“ 

0.0034 

0.91 1 5 

3 

0“ 

0.7548 

0.834J 

BaSO, 

16 

3 2“ 

0.7581 

0.8357 

Dust"* 

60“ 

0.8335 

0.8713 

0“ 

0.9.174 

0.9438 

64 

12“ 

0.9185 

0.9444 

60“ 

0.9437 

0.9563 

0“ 

0.9487 

0.9595 

256 

12“ 

0.9493 

0.9599 

60“ 

0.9648 

0.9685 

I 
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C.  Discussion  of  Results 

1)  relation  lietv/een  Imaginary  index  and  absorption  coefficient. 

Tile  discii.sion  In  appendix  D sliows  that,  in  general,  the  plausible  hypo- 

V 

thetical  relation  a Nv(4irK/X)  is  badly  violated.  Tiie  results  in  Tables 
2 and  7 reinforce  this  conclusion;  the  values  of  a'  vs.  a show  discrepancies 
of  a factor  of  two  or  three.  For  example,  following  the  reasoning  in  the 
latter  part  of  appendix  D,  one  could  attempt  to  predict  a value  for  k of 
carbon,  from  the  correct  calculated  v;ilue  of  a given  in  table  2,  for 
A = 0.6  pm.,  using  tlie  plausible  relation  above.  From  table  2,  sample  2, 
for  carbon,  v = 3.142  x lO"^  (pm)^,  N = 3.744  x lo"^  (pm)“'^,  a = 5.205  x lO"'"* 
(pm)  The  hypotlietical  relation  above  yields  k = A a/4ir  Nv  = 0.21,  rather 

than  tlie  correct  (input)  value  k = 0.6. 

2)  Relation  between  (K,S)  and  (K,S) . 

Tables  3 and  8 display  tlie  values  of  K,S  calculated  from  F,qs.  (152), 
using  the  accurate  multiple  scattering  values  of  (R^,R^).  Kxperlinental  ists 

P 

often  hypothesize  that  these  values  should  be  equal  to  (K  “ 2oi,  S t 20  ) , 
regardless  of  the  beam  angle  of  incidence  or  of  the  optical  depth  of  the 
sample.  The  calculated  results  in  these  tables  show  that  (K,S)  depend 
somewhat  on  both  0^  and  d.  However,  In  the  cases  of  interest,  samples  2 
and  3,  setting  2a  equal  to  any  of  the  K values  leads  to  an  error  of  less 
than  20%,  which  is  probably  acceptable  at  present.  But  the  values  of  S 
seem  much  closer  to  P_  itself  rather  than  to  20  , the  expected  correlation. 

The  only  large  discrepancy  between  K and  2it  occurred  at  both  wavelengths 
for  the  pure  BaSO^^  sample.  I’ractically  speaking,  this  is  irrelevant,  since 
the  absorption  due  to  HaSO^^  In  a mixture  is  probably  negligible  in  most 
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actual  experimental  cast!S,  and  is  certainly  negllRible  in  the  model  samples 
treated  here. 

The  computer  c.-ilculatlons  are  accurate  enoiiidi  so  that  the  above-mentioned 
discrepancies  (S  - p_  rather  than  2p_,  K >>  2ct  for  the  very  weak  absorber 
BaSO^^)  are  probably  rcsal,  not  just  artifacts  of  the  arithmetic.  If  so, 
it  may  be  conjectured  that,  since  the  model  samples  all  contain  closely 
packed  particles,  which  are  strong  reflectors,  the  origin  of  these 

discrepancies  may  just  be  what  is  discussed  in  appendixes  K and  I,  namely, 
the  failure  in  principle  of  the  standard  radiative  transfer  equations  for 
closely  packed  strongly  reflecting  or  absorbing  media. 

3)  Near  field  vs.  far  field  phase  functions. 

The  underlying  theory  for  the  near  field  phase  function  is  discussed 
in  appendix  A.  T hies  5 and  6 present  comparisons.  Although  the  angular 
dependence  of  the  near  field  phase  functions  is  significantly  different 
from  tliat  of  tlie  corresponding  far  field  ones,  the  total  reflectances  were 
not  very  different,  for  the  samples  considered.  Table  2 shov/s  that  the 
packing  fraction,  the  ratio  of  tlie  ac;tual  volume  of  an  average  particle  to 
the  average  volume  subtended  per  particle  in  the  mixture,  equal  to  Nv,  was 
of  the  order  0.')  for  these  samples.  It  is  possible  that  larger  packing 
fractions  would  produce  greater  difference.s  between  near  and  far  field 
pr(>dictions  of  reflectances;  this  remains  for  future  work. 

4)  Reiclimaa  tlieory  vs.  full  multiple  scattering  tlieory. 

2 

Tlie  Kelciim.in  tlieory  is  discussed  in  appendix  C.  Tables  4,  6,  and  9 
give  comp.arisons . 'fho  discrepancies  are  generally  around  10%,  so  this  two- 
flux  theory,  which  takes  into  account  the  incident  beam  direction,  could 
in  principle  be  used  to  evaluate  (K  = 2cx,  S = 2(1  ) from  reflectance  data. 
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'llie  formiilafi  aro  i.  i'  complicated  than  fhose,  of  the;  KM  theory,  and  more 
data  i ;;  nf!odcd  for  e.ach  attiuiipted  evaluatLon,  Imt  mort:  information  about 
Mie  theory  parameters  is  obtained.  back  inf;  an  error  analyr.is,  it  Is 
difficult  to  estimate  whei.'uer  the  Ucichman  tlieory  v;ould  be  more  accurate 
in  iiract  Icab  predictions  of  (K,S)  than  the  KM  tiieory. 

IV.  Tlieorctical  Studies 

Tlie  fundamental  theoretical  studies  which  were  done  during  this 
investigation  centered  on  the  problem  of  multiple  scat: ‘ring  in  closely 
packed  strongly  absorbing  or  strongly  reflecting  media.  Two  major  questions 
arise  in  tiiis  problem.  Tlie  first  question  asks  whether  the  standard  differ- 
ential equation  of  r;idiative  transfer  should  be  applicalilc;  to  sucli  media. 

The  second,  wliich  arises  only  if  the  first  is  aiiswered  negativt-ly,  concerns 
wlietiier  <a  suitable  equation  of  transfer  can  be  found  and  solved. 

In  this  work,  tlie  first  question  was  answered  negatively,  and  the 
second  affirmatively.  'J'he  detailed  derivations  and  calculations  are  pre- 
sented in  a|)pendixes  K,  A,  and  I.  In  app>-‘ndix  1,  it  is  shown  that  tlie  new 
radiative  transfer  equation  .sliould  yield  sij>n  i f i c.antly  different  results 
than  the  standard  equation,  in  many  cases  of  practical  interest. 

Since  tlie  proposed  new  equation  of  radiative  transfer  may  have  practical 
consequences,  it  is  Important  to  achieve  a thoroug.h  understanding  of  the 
logical  steps  which  lead  to  the  equation.  I'or  the  interested  reader,  the 
logic  patliway  starts  witli  aijpenilix  K,  then  j;oes  to  appendix  A,  then  finally 
to  appendix  i. 

In  appendix  K,  a model  one-dimensional  random  scattering  medium  is 
constructed.  flie  medium  consists  of  a stretched  fric.tiotiless  string,  with 
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I 

I 

I 

I 

I 

I 


randomly  but  homogt'.noiisly  placed  Insertions  of  segments  of  a frictional 
string,.  Tliis  medium  contains  analogies  to  every  major  feature  of  a 
3-dimensional  random  scattering  mediirm.  The  wave  motions  are  of  course 
the  transverse  oscillations  of  the  string.  Tlie  exact  equations  for  the 
intensity  transfer  tlirough  the  medium,  which  are  the  complete  analogs  of 
the  equations  of  radiative  transfer  in  a 3-dimensional  medium,  are  derived 
and  solved  without  approximation.  Tlie  equations  turn  out  to  be  difference 
rather  than  differential  equations;  they  reduce  to  differential  equations 
of  standard  (two-flux  theory)  form  only  if  the  Inserted  scattering  segments 
are  very  weak  absorbers  and  reflectors.  Tlie  exact  solution  of  the  correct 
difference  equations  for  either  strongly  fibsorbing  or  strongly  reflecting 
Segments  is  quite  different  from  the  incorrect  differential  equation 
solution,  and  the  expressions  for  transmittance  and  reflect.ance  are  also 
quite  different. 

In  appendix  A,  the  derivation  of  tlie  usual  3-dimensional  differential 
equation  of  radiative  transfer^  is  done  in  an  unusual  manner  which  is 
partly  analogous  to  the  derivation  done  in  appendix  K for  the  one-dimensional 
case.  Indeed,  a difference  equation  (Eq.  A6)  arises  in  tliis  derivation  for 
the  3-dimensional  case,  and  it  is  necessary  to  require  yP,  <<  1 in  order  to 
achieve  the  usual  differential  equation  (A9) . Here,  y is  the  extinction 
coefficient  of  tlie  medium,  and  P,  is  the  average  interparticle  separation. 

If  yP  <<  1,  then  the  change  in  intensity  acrtiss  a mcmoparticle  layer  is 
negligible.  But  if  yP  is  not  very  small,  there  will  be  enough  change 
across  a monopart  icle  layer  that  the  difference  equation  results  '.'ill 
disagree  with  the  differential  equation  results.  Tlie  value  of  yP  which 
should  be  used  as  a break  point  between  the  two  equations  depends  of  course 
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on  the  accuracy  desired  in  a calculation.  If  accuracy  1%  is  dor.  i red, 
then  the  difference  equation,  rather  than  the  standard  differential 
equation,  should  probably  be  used  if  yt  > 0.01. 

Although  a near  field  rather  than  a far  field  phase  function  should 
be  used,  in  principle,  for  closelj'  packed  media,  the  computational  results 
quoted  in  section  TV  imply  that,  practically  speaking,  the  two  are 
equivalent . 

In  appendix  I,  the  difference  equation  analogs  of  Eqa.  (A18),  (A20) 
are  derived.  These  difference  equations  are  applicable  to  problems  with 
(x-y)  translational  invariance  in  which  azimuthal  dependence  is  irrelevant, 
such  as  the  problems  which  were  treated  in  this  work.  Methods  of  solution 
are  discussed;  a combination  of  the  many-flux  method^'^  and  the  doubling 
method  is  suggested.  It  is  also  conjectured  that  present  Monte  Carlo 
methods  might  already  be  solving  the  difference  equation,  since  these 
methods  treat  the  radiative  transfer  in  discrete  steps,  collision  by 
col  1 ision. 

Also  in  aijpendix  I,  the  two-flux  difference  equations,  analogous  to 
the  KM  equations,  arc  developed  and  solved.  TJie  results  are  contrasted 
witli  the  corresponding  KM  theory  results.  By  means  of  two  realistic  examples, 
it  is  shown  tliat  the  discrepancies  between  the  two  theories  probably  have 
practical  significance.  In  fact,  one  of  tlie  examples  involves  the  realistic 
BaSO^^  powder  model  which  was  used  in  the  computations  in  this  work,  llie 
parameter  values  in  table  3 lead  to  P_S.  = 0.283,  Indicating  that  the 
difference  equations,  rather  tlian  the  differential  equations,  probably 
should  have  been  used  in  all  tlie  computations  which  were  reported!  Time 
did  not  permit  thin;  it  is  left  for  future  investigation.  However,  the 


principal  source  of  terror,  in  obtaining  values  of  (K  2a,  S 3 23  ) from 
difliise  reflect  anc<*  data  by  t.ieans  of  the  KM  theory  relations  (B2)  , lies 
in  the  neglect  of  the  incident  beam  an'/,le  dependence;  this  nej'lect  v/ould 
not  In;  removed  merely  by  use  of  the  difference  equation  analog  to  the  KM 
equations.  And,  of  course,  multiple  scattering  has  nothing  to  do  with  the 

7 

large  error  in  the  hypothetical  rf:lation  a = Kv(4ttk/A)  which  has  been  used 
to  obtain  values  of  k from  the  value  of  K = 2a  determined  as  above. 
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APPENDIX  A 

RADIATIVE  TRANSFER  THEORY  AND  NEAR  FIELD  PHASE  FUNCTIONS 

Im.'ip.ino  n beam  of  radiation  incident  at  aone  anj;le  0 with  respect 

o 

to  the  normal  to  the  plane  surfaett  of  a homo}>eneoiis  random  scattering 
medium.  If  the  beam  has  a cross-section  broad  enough  to  subtend  a very 
large  n»imber  of  particles  in  the  medium,  then,  for  purposes  oi  analysis 
the  actual  collection  of  particles  of  different  optical  types  and  sizes 
may  be  replaced  by  a collection  of  identical  "average"  particles,  as 
follows . 

Define:  N = total  number  density;  = number  density  of  type  i 

particles;  f^(r)dr  = fraction  of  type  i particles  with  radii  in  (r,  r+ 
dr),  jdr  f^(r)  = 1;  o^(r,s,s')  = differential  scattering  cross-sec ti<jn 
of  a particle  of  radius  r,  type  i,  from  direction  s'  to  s;  > 

o^  (r) > (r)  S (extinction,  scattering,  absorption)  cr<'  ;-sections 

SC  cl  Hi)S 

of  such  a particle.  Then  the  average  particle  has  the  following  differ 
ential,  extinction,  scattering,  and  absorption  cross-sections: 


o(s,s') 


dr  f^  (r)o^.  (r,s,s' ) 


ext ,sca ,abs 


N-‘  I N, I 

H N ^)^N.jdr  f.(r)o\r) 
i •’ 


(Al) 


ext , sea , ahs . 


Here  and  in  what  follows,  s,  s'  are  unit  vectors.  Imagine  the.  medium 
divided  into  identical  cubical  cells,  wltli  side  length  equal  to  the 
average  interparticle  spacing  f.,  witii  = 1/N.  Each  cell  contains  one 
average  particle,  on  the  .average,  located  randomly  wltliin  the  coll.  As 
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lonf,  as  the  randomness  in  particle  location  extends  over  a wavelenj'.th  or 
more,  radiation  intensities  may  he  used,  rather  than  phase-sensitive 
amplitudes . 

Consider  the  cell  with  center  at  position  x,  and  radiation  on  tliis 

cell  in  direction  s in  dfl  . The  orientation  of  a cell  is  immaterial,  so 
~ s 

it  may  be  chosen  with  two  of  its  faces  perpendicular  to  s,  for  any  s. 
Therefore,  the  incident  power  on  the  cell  is 

P.  df>  = f.2  I(x-^s/2,s)dSl  , (A2) 

in  s ~ ' s 

where  I(x-i.:5/2)  is  the  intensity  on  the  incident  cell  face,  which  his  a 
noimal  parallel  to  s,  and  a center  at  x - is/2.  This  beam  encounters 
the  average  particle  somewhere  in  the  cell,  and  scatters  anti  absorbs 
radiation.  The  power  lost  out  of  dfl^  is  thus 


P dT  = o T (x-ts/2)d:>  . 
ext  s ext  ' B 


(Al) 


On  the  other  hand,  part 
directions  gets  scattere 

P dP  = dn  [ 


of 

d 

dP 


sea  s s 


the 

into 


power  from  incident 
dP^  by  the  particle 

,s’)I(x-Hs'/2,s') 


beans 
; this 


in  all  other 
ix'wer  is 

(AA) 


wlicre  l(x-is'/7)  is  the  intensity  at  an  incident  cell  face,  with 
center  at  x-f.s'/2,  normal  parallel  to  s'.  If  designates  the 

power  emerging  from  the  cell  in  through  the  exit  fiice  centered  at 

X + 7.s/2,  then 

P dP  = I(xbfs/2,s)dP  . 
out  S ' - - s 


(AS) 
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P^nerw  halniice  requires  P = P.  ~ P +1’  .or 

out  in  ext  sc;i 


I(x+J,r./2,s)  = 1 (x-J-s/2)  ! Y^. (Alt) 


-1 


, P (k  , <j  ' ) I (x-7.<j  '/2  , s ’ ) , (Afj) 


wlierc 


Y = Na  . = extlm  tion  coefficient 
ext 


p(s,s')  = (411)0  (s  ,s  ' ) /a  ^ pliase  function. 
~ ext 


(A7) 


(A8) 


, -1 

Remember,  7,  = N . Now  make*  a Taylor  expansion  at)out  the  point  x.  If 

the  chanp.e  in  I over  distance  £/2  is  small,  (which  will  be  tlie  case  if 

Y^  <<  1),  then  terms  proportional  to  Y^^'IY  he  neftlected,  and  Kq. 

r, 

A6  reduces  to  the  standard  radiative  transfer  (K'l')  equation  , 


y *s-V  l(x,s)  = -I(x,s)  + (4it)  ' dSl  , p(s  ,s' )I  (x,s' ) 


(A9) 


Tlie  quantities  (e.,f.),  tiie  (absorption,  sentterir,!’)  coefficients,  arc- 
defineii  by 

a B = No^^^.^.  (AlO) 

Note  that  Y = a+B,  and  that  the  phase  fund,  ion  is  normalized  so  that 


(An) 


dP^  p(s.s’)  ^ = B/y  ^ . 


(All) 


the  (sinj'.le  scat  tc-r  i up,)  alhedo. 

This  metliod  of  derivation  makes  it  clc.-ar  just  liow  the  scalterinf; 
properties  of  tlie  individual  particles  in  a random  mixture  enter  into 
tlie  radiative  transfer  equation. 

An  interesting  point  arises  for  close-packed  media,  where  a cell 


JT‘  • ;r  '■ 
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is  not  imich  bigger  than  tlie  aver.’gc  scattering  particL;'.  Tlie  above 

method  of  derivation  makes  it  cl  t that  the  sca*'toror  in  the  next  cell 

(in  direction  s)  seta  the  near  zone  and/or  induction  zone  fields, 

rather  than  the  radiation  zone  fields,  from  the  scatterer  in  tlie  given 

cell.  I’andomr.e^s  cf  scatterer  location  extending  over  distances  > X 

still  allows  use  of  intensities  rather  than  phase-sensitive  quantities, 

but  o(s,s')  should  be  replaced  by  o (s,s*)  in  the  expression  f(>r 

P , where  a (s,s*)  is  a suitable  near-field  cross-section, 
sea’  near  ~ ~ 

In  this,  work,  o (s,s')  was  defined  as  follows.  Consider  tl... 
near  ~ - 

Mie  scattering  problem,  with  a spherical  scatterer  centered  at  the 
origin  of  coordinates,  and  a plane  electromagnetic  wave  of  unit  power/unit 
area  incident  in  direction  s'.  Let  S(x,x,s')  ^ ti  i.:e-averaged  Poynting 
vector  of  the  scattered  fields  at  distance  x in  direction  x.  Tlien 

S(x,P.,s')  = E (x,t)x  h (x,t). 

~ Att  'scA'  ~sca  ~ 


The  near  field  scattering  cross-section  is  defined  in  this  werl;  by 


2 

o (s,s')dfi  H X s • S(x,  s ,s' )dP.  . 
near  ~~  s s 


(A12) 


Note  that  as  x •>  a (s,s')  •>■  o(s,s'),  tbe-  usual  d i f feri'nt  ial  cross- 

near  ~ ' ' ' 

section  defined  in  terms  of  the  .asymptotic  fields.  The  distance  x was 
taken  to  be  equal  to  some  factor  times  the  particle  radius  r,  for  each 
particle  in  a mixture,  in  this  work.  Explicitly,  defined  by 


a 

near 


N 


-1 


I N, 


dr  f^(r)o^ (Fr ;r ,s  ,s' ) ; 


(All) 


J 
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that  is,  the  ficlcis  u.sed  In  i-valiiatin^  tlie  cross-scr t ion  of  a particle 

of  radius  r .m-c  those  at  a distance  J'r  from  the  particle's  center,  since 

these  are  approximately  the  ones  incident  on  the  nei  p.hhor in};  particles. 

If  N = l/l^  is  the  overall  number  density,  v H N ^Yn.v.  is  the  volume  of 

. 1 1 

an  average  particle,  then  the  factor  F should  be  chosen  as 

F = = (1/Nv)’'^^  . (A14) 

Note  that  F > 1 as  long  as  there  are  air  spaces  between  particles,  v;hich 
there  must  be  if  the  medium  is  to  be  a random  scattering  medium. 

It  is  straightforward  but  nontr Iv'ial^^  to  show  that  the  above  defini- 
tion of  the  near  field  cross-section  yields 


a (s,s’)dSJ  = o , 
near  - ~ s sea 


(A15) 


4ti 


where  o is  the  conventional  asymptotic  total  scattering  cross-section, 
sea  ' 

independent  of  x.  Since  any  absorption  by  the  scatterer  is  strictly 


localized  to  regions  within  the  scatti'rec,  conservation  of  enerjjy  shows 

that  a + o , = o » also  independent  of  x.  Define  the 

ext  sea  abs  ext  ‘ 

near-field  phase  function  by 


P (s,s')  4ll  No  (s,s')/y,  (4ti) 
near  - ~ neai  - - 


-1 


Ait 


d,Q  p (s,s')  = m . 
s near  ~ ~ o 

(A16) 


Tills  phase  function  has  a different  angular  depemlence  than  the  far 
field  one,  although  its  normalization  is  the  same.  An  appropriate 
radiative  transfer  equation  for  close-packed  random  scattering  media  is 


then 
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Y ^ s-VT(x.s)  = -I(x,r)  + (4ti)  > § ' ) I (x,  s ' ) , (A17) 


as  lonf;  ns  y9.  <<  1.  If  yH  < 1,  a simple  dlf fercMitial  t‘.c]uation  c.'nnot 
be  obtained  from  F.q.  (A6)  ; this  case  is  considered  in  appendices  E and 


For  problerns  with  (xy)-LransIatlon  invariance,  which  are  the  only 
oi'.es  treated  in  this  work,  Eqs.  A9  and  A17  reduce  to  the  forin 


± XTI 

P = “YT  (?.,]!,(!))  +*^1  ) * (A18) 


-1  o 


where  (>=cos  p,(}>)  are  the  polar  anp.les  of  the  unit  vector  s with  respect  to 
the  z-direction.  Since  only  total  fluxes  (transmitted  and  reflected) 
are  of  interest  ia  this  work,  the  relevant  radiative  transfer  equation 
is  obtained  by  integrating,  Eq.  (A18)  over  all  azimuth  angles  (} . 

He fine 


T(z,ii)  r (2ir)  d(>  I(:-,p,(J));  p^‘^’^(y,p')  = ^ d<|i  p (p  ,<> , P ' , v ' ) 


(A19) 


Tlicn  Eq.  (A18)  reduces  to 


P = -yT(z,p)  + -p' j dp'  p^”\p,p*  )T  (z,p' ) 

-1 


(A20) 


which  is  the  standard  azinuth-indcqiendent  radiative  transfer  (equation 
for  problems  with  (xy)-translat  ional  invariance.  Hero,  p^°\p,p')  is 
formed  from  cither  the  ne.ar  field  or  far  field  phase  functions  discussed 


above . 
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If  thcro  is  a uniforni  incident  bean,  siicb  that  I (O.p.A)  = 

j II 

7tF6(vi-jj  )6(<>-<}>  ),  tlu-n  T.  (0,|()  ^ ),  and  T.  (k,ij)  = 

o o in  / o in 

1 — Y Z / Ll 

- l’6(vi-p  )c-  ' . Subst  itutiiif'  this  into  (A20) , with  l(z,p)  = I.  (z,p)  + 

^ O 11 

l^,i(z,p),  yields; 


ai 


V - 


di 


-Yldi(^-.h)  + 2 Y 


1 

r 

dp'  p 
-1 


(o) 


(iMi'  )Idi  (7-,p'  ) + F 


(O),  , 

P (P,P  )e  , 


(A2i) 


which  is  the  relevant  radiative  transfer  equation  for  the  diffuse  inten- 
sity I^^(z,p).  It  is  thin  equation  which  v;as  solved  accurately,  usinp, 
computer  alp.orithms  based  on  the  doublinj>  metliod  of  Van  de  Hulst  and 
Hansen,  for  media  of  various  thicknesses  d (optical  thickness  t = y’)- 
Both  near  and  far  field  phase  functions  were  used  and  compared. 
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API'EKDIX  P. 

KUBKl.KA-;nP«  THHOKY 

The  Kubp] k.'i-Munk  (KM)  theory^  is  equivalent  to  the  Schuster  tv/o- 
12 

flux  theory  of  rc.cllation  transfer  In  plane  parallel  hoino{;eneous  random 
sentijrinp,  media.  Tn  Eq.  (A20) , assume  that 


= l^(z),  II  > 0;  T(z,p)  = I_(z),  p < 0 

and  Intof'rate  Eq.  (A20)  over  0<vi<l,  and  over  0>p>  -1.  This  yiel<ls  the 
KM  equations 


d4> 

dl  --(K+S)(.^+S! 


dz 


- (K+S)1>  - S«. 


1 1 


where  K 


2n,  S - 2P_;  H ^ 


clp 


dp'  p^°\p 


,iTi');  S H j 


(VA) 


1 

dp  p T_^  (z)  = 


I , - (magnitude  of)  fluxes;  and  ^ y = a+B . + P . Note  ?hat  S is  twice 

the  hack-scattering  coefficient,  and  K is  twice  the  absorption  coefficient, 

13 

just  as  in  the  usual  ad  hoc  derivation  of  the  KM  equations  . In  this 

work,  the  sample  liad  thickness  d,  and  was  placed  on  a black  background  of 

reflectance  0.  The  solution  of  Eqs.  (PI),  subject  to  the  boundary 

conditions  ('^  / 'i'  ) = R,,  ('5  / i' . ) , = E = 0.  yields  the  KM  relations 

+ z=o  d - + z~d  g 


K 

S 


(1-R  ) R 

— , .Pd  = p £n 

^ “ 1-R^ 


/l-R.R 

Cl 

iT^Tir 

* d « 


(P2) 


whore  R^  is  the  reflectance  of  an  lnfinit(*ly  thick  san)]>lc.  Tliese  relations 
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have  often  been  iiscil  Lo  evn]>iate  K,  S,  fxom  measureil  values  of  R,,  R . 

cl  o> 

In  the  usual  ad  hoc  derivation  of  tlie  ihM  tlieory^'\  llie  Inciclcuit 

illunination  is  supposed  to  bo  eitlicr  diffuse  and  isotri'pic,  or  a beam 

at  f = 60"  (p  = l/:i). 
o o 

Koto  v.'cll  that  the  two-flux  approxiirat  ion  to  the  full  RT  Kq.  (A20)  is 
very  crude,  and  that  all  directional  effects  which  would  be  caused  by 
bear,  incidence  at  an  arbitrary  anp.lc  0^,  rather  than  by  diffuse  isotropic 


incidence,  are  lost. 


KF.tCHM/V','  'lIii:OFY 


2 

This  two-flux  theory,  dc/rl oped  hy  Kcichiiar)  > stnrLs  from  Kq.  (A21) 
rnther  than  (A20) , so  it  represents  an  atteri|)t  to  include  incident  beam 
directional  effects.  In  Kq.  (A21),  assume  that 


Idj(j^.l')  = h>0;  = T_(^).  P<0;  K ^ 1; 


and  intC'j’.rate  Kq.  (AlO)  over  0<p<l , and  over  0>|i>  -1.  Tit  is  yields 


d«>  -Y7./II 

— - = -(K+S)<f.  + S?-  + yP.(p  )c  " 

d z + ~ + c) 

(Cl) 

d<  -Yx/p 

--  ■ = (K+K)0  - S4'  - yT>  (h  )e  ” 

dx  - + - o 


where  K 2a,  S H 2B_,  ~ 1^/2,  .just  as  in  Appendix  li;  and 

1 


1\  (P  ) " 7 
+ o A 


f'di  P^)- 


o 


(C?) 


Kquatif'iin  (Cl)  are  very  similar  to  the  KM  equations  (111),  except  that 
here,  , are  diffuse  fluxes  only;  the  incident  flux  is  ~ 

^ -vz/p  p^  -Yx/p^^ 

dp  p --  6(p-p  )e  = — r e . Boundary  rondit  ions  are  0 (d)  = 0 

J 2 o / 

o 

(l)lack  haclj’.rounil)  , 0|(0)  = 0.  The  reflertant  e is  (f  (0) (D) 

for  a sample  of  thickiuuts  d. 

In  p.eneral,  the  tiolution  of  these  equitions  with  ther.e  boundary 


conditions  yit'lds 
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o 


2 2 

where  L H K + 2KS, 


SD^sinh(Ld)  + l.h  exp-ycl/M^^ 

" ’^’sf  'co7ir(rd)~l^“l7  sinlT(i,dT“ 

and  D H y[(K+S  ± y/h  )P,  + S)’  ] / (I,^-y . 

~ 0+7  O 


(C3) 


In  the  special  case  that  y (‘T-gle  of  incidence  - 60°)  and 

B,  = B =4-,P  =P  = io/6,  (isotropic  scatterers  will  yield  this), 

Eqs.  (C3)  yield  just  the  KM  relations  (B2)  for  K/S  and  Sd  in  terr.s  of 
and  R^.  In  all  other  cases,  the  expressions  for  K/S  and  Sd  are 
different . 

This  means  that  the  original  KM  theory  should  be  just  as  good  an 
approximation  as  this  modified  theory  if  0^  = 60°,  and  if  the  scattering 
is  isotropic;  but  that  the  KM  theory  should  be  Ics  : accurate  than  this 
modified  theory  for  beam  Incidence  at  other  angles,  and/or  for  non- 
isotroplc  phase  functions. 

In  this  work,  the  predictions  of  R^,  R^  from  tills  raodifie’d  theory 

wave  checked  against  those  from  the  accurate  doubling  method  solutions 

of  the  RT  equations,  and  were  fiuind  to  agree  within  157  for  all  samples 

? 

which  were  modelled.  Similar  ag.rt'ement  bar.  been  noted  by  Reiehman  . 
Equations  (C3)  may  be  rewritten  in  the  follouiug,  form,  with  cxp(-yd/p^)  : 0 


8R  (2R  la,  - a ) 
K_  _Y.  J_  _ 

(1+R  )^-4p‘'(l-R 

K OK 


(f:4) 


-2I,d,  ...  ^ 1 -2Ld, 

e.  )/(l  + I<  e ), 


where  K/S  (1-R  ) /2R  ; a^  ^ P^(p  )[(l  ( K/S)p  ' (1  + K/2S)  1 > )i  P (p  ), 

K K 7.  + O O O - C)  ' 

r = (1 

For  samples  of  optical  thicknes.s  yd  > 4,  this  is  a very  g.ood  approxima- 
tion. In  principle,  measurement  at  ;i  g.iven  p of  R and  R,  at  three 

O a*  Cl 


i 


37 


I 

I other  d,  will  allow  Inversiiou  of  Eq.  (C4)  to  obtain  values  of  the  four 

unknowns  , L,  a_.  In  turn,  knowle ’“e  of  these  four  quantities 

I allows  determination  of  the  Mie  theory  parameters  (e.,P^,y),  or,  since 

N is  known,  the  cross-sections  a . , a , o ^ for  the  averaee  particl;'. 

P in  the  medium.  If  the  particle  types,  number  densities,  and  size 

II  distributions  are  known  reasonably  well,  Mie  tiieory  computer  aljjorlthms 
could  be  used  to  zero  in  on  values  of  the  real  and  imaj;inary  refractive 

H indices  of  each  type  of  particle  in  the  medium.  This  possibility  was 

not  pursued  in  this  work,  but  it  seems  to  merit  further  study. 
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AIM')  : 'IX  U 

RK1.AT10N  BKl  i.'Ki'N  k AND  « 

Con3ick*r  a randon  scattcTing  madiiim  compospd  of  N part  iclos/iini  f 
volume.  For  convenience  in  what  follows,  let  all  th<-  particles  liave 
the  same  radius,  volume  v = 4nr^/3,  and  the  same  refractive  inde/-:, 
m = n-iK;  the  discussio t presented  below  can  be  extended  trivially  to 
Include  the  general  case. 

As  long  as  the  arrangement  of  particles  is  homogeneous,  )>ut  the 
randomness  in  the  location  of  each  particle  is  greater  tlian  a wavel  csij’.th , 
X,  the  condition  for  applicability  of  the  standard  radiative  transfer 
equations  (A9)  are  met,  and  the  absorption  coefficient,  a,  of  the  medium 
is  defined  in  terms  of  the  Mie  theory  absorption  cross-sect  ion , n , by 

cli'S 

a “ No  , . (Dl ) 

abs 

The  derivation  of  the  KM  equations  in  Appendix  l’>,  and  the  Ueicbman 
equations,  in  Appendix  C,  makes  it  clear  that  the  KM  absorption 
coefficient  K is  defined  by  K = 2u. 

How  imagine  that  this  collection  of  i)articles  is  fused,  with  the 
particles  distorted  as  necessary  and  assembled  so  that  tlu're  are  no 
longer  any  air  spaces  betv'ern  thc-ra.  Then  this  collc'ction  of  particles 
is  indeed  a uniform  material  medium,  possessing  as  a whole-  a rc-fractive 
Index  m = n - ix.  llie  absorption  coefficient  = extinction  coefficient 
of  this  medium  is  obviously 

a = Atix/X,  (D2) 

m 


.■2T ’.r-** 
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This  follov/s  from  the  behavior  of  a plane  v..ive  exp  i(o)t-2Tim  z/X) 
travelinp.  In  any  diroctioa,  say,  tlic  z-d  i i C’Ction . 

In  some  work,  the  plausible  Iiypothesis  has  been  riade  that  u should 

in 

be  related  to  a simply  by  the  ratio  of  the  total  volume  occupied  by  the 
random  scatter  inf;  nuidiurn  to  that  occupied  by  the  same  total  number  of 
particles  arranged  as  a imiforn  material  mod Lum. ^ ^ That  is,  the 
hypothesis  is 

? 

a = Nv  = Nv(4inc/X).  (D3) 

That  this  hypothesis  is  plausible  follows  from  the  reasoning  that 
the  amount  of  energy  absorbed  by  a particle  sihould  be  proportional  to 
Its  volume,  for  a given  incident  radiation  field.  However,  it.  is  clcar 
that  the  absorption  by  a volume  element  of  a substance  is  proportional 
to  the  local  Intensity  inside  that  volume  element,  not  to  the  incident 
intensity  on  the  element.  If  the  (average)  particles  in  a random 
scattering  medium  arc  compressed  so  tightly  that  the  medium  becomes 
a true  ma'  rial  medium  (lU'  boundaries,  no  diffuse  scattering),  then 
the  local  Intensity  in  each  vt'lum.e  element  which  was  subtended  by  a 
g.iven  imrticle  is  completely  iliffc^rent  from  the  local  inti‘n,;ity  which 
existed  inr.idj'  the  ( orrt-sponding,  volur'e  element  in  the  random  scattering 
medium.  In  g.oinj’,  from  r.catter  ing,  medium  to  material  medium,  local 
Intensity  is  not  invar iint,  which  it  would  have  to  be  in  order  that 
the  hypothesis  be  valid.  Note  that  it  is  the  single-scattering  local 
intensity  that  is  considered  here,  because  the  absorption  coefficient 


/.o 


of  a random  scaM  cTinj;  mod  luni  is  defined  in  terns  of  a layer  of 
partic,1-'S  so  thin  that  multipic  scattering  in  that  layer  may  he  nt'p.l  tfct  ec' . 
Tlierefore,  the  quostiyp  of  whether  the  local  intensity  chanpes  can  he 
answered  by  single  scattering  theory,  e.g. , Mie  theory  for  the  sptierical 
particles  considered  in  this  work. 

Combining  Eq.  (01)  and  the  bypotbesis,  Eq.  (Oil),  yitlds  tbe 
equivalent  bypothet ical  relation 

7 

o , = Atiitv/X.  (OA) 

cl  i) 

This  hypothetical  relation  is  expressed  entirely  in  terms  f)f  Mie 
theory  parameters,  and  nay  be  chi-cked  by  that  (lieorv  alone.  .Several 
examples  will  serve  to  show  that  the  byiiotliesis  is  valid  only  for  small 
particle.';  whose  refractive  index  i.s  close  to  1,  and  badly  violated  for 
other  cases. 

lA 

Van  de  Ilnlst  gives  the  following  approximate  expression  for  a , , 

abs 

valid  for  (?irr/X)  ^<<1  , for  any  in,  for  a particle  of  radius  r: 


(1)5) 
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case  1)  (2TTr/A)^  <<  1 , ii  = 1 , < <<  1 . Tht>n  only  the  first  term 

7 2 

in  Eq.  (D'>)  cont  r ibiit  i!S ; Im((n  +2)  = -7.yJT>,  aiv’  (t  = 

at).'! 

3 

(4Tir  /3)(4ti/A)k  - 4ttkv/A,  in  n>', reorient  with  the  iiypot  hesi  s , Eq . (D4). 

case  li)  (2iir/A)^  < 1 , n = 2 , if  <<  1 . Then  Eq . (1)5)  yields, 
after  some  ar  iLhiretie  , 

i 

'’ahs  = (^•’’<v/A)(~  + — (2TTr/A)2).  (l)f>) 

It  is  seen  that.  If  (2-nr/A)  <<  1,  the  hyiiothesis,  Eq.  (1)4),  is  violated 
by  a factor  of  2;  that  is,  usinj’,  Eq.  (D4)  and  an  experimentally  deter- 
mined = a/N,  one  would  predict  a value  of  r only  half  an  lar<>e  as 

the  actual  value.  I’or  (2-nr/A)  x 1,  Eq . (I)’-)  seems  to  ap.roe  with  the 
hypothesis,  but  tin-  iiep.l  -cted  terms  of  0(2nr/A)''  and  hip.her  in  Eq . (1)5) 
make  Eq.  (I)f>)  Invalid  for  (2rr/A)  : 1. 

The  jjeneral  case,  V7ith  both  2TTt/A  and  m allowed  to  have  any  reason- 
able values,  cannot  be  analyzed  except  by  use  of  the  full  Kie  tlieory.  for 
example,  in  the  limiting  case  of  very  lar;’i’  .scatterers  (2nr/A  >>  1), 
the  asymptotic  relation  from  Nie  tluory  is 

o , = 7tr’(l-m  ) - (1-m  ),  0<ui  ^1.  (D7) 

abs  o 4r  o ~ o“ 

The  hypothesis,  Eq.  (1)4),  would  thus  reijuire 

]f))'i73A  - (l-w  ), 
o 

which  is  clearly  very  badly  violated  for  2nr/A  >>  1. 

Standard  Mie  theory  computer  alp.orithms,  such  as  those  used  in  this 

work,  can  accurately  and  easily  preilict  values  for  o , , o , o , ..., 

^ ' abs  sea  ext 


I 
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at  any  A,  p.lven  i)  tlu;  Dptical  ccnritants  t'>‘‘  tyP'  i pari  id 

in  tiio  nccliun  to  be.  modi’ll  d,  and  ii)  tlio  .‘-.izc  d isti  Ibnt  ion  f^(r)  of 
each  type’  of  particle.  I'n  fo?  t unatoly , it  is  not  so  easy  In  invert  these 
calculations,  not  even  in  the  r.iii’plest  case  of  a inonod  i ripein;e  collec- 
tion of  one  optical  typo  of  particle.  In  tliis  simplest  case,  (or  ex.imi'lo, 
suppose  the  real  refractive  index,  n,  i r.  Inown,  and  the  radius  r is 
known,  for  a constituent  particle.  If  tlieri-  is  at  least  an  order-of- 
maf;nitndo  feel  inp.  for  v;hat  tlie  value  of  the  imap, inary  index  < oup.hl  to 
be,  then  the  standard  computer  al por i tlims  can  he  used  to  search  around 
the  expected  value  of  k until  tluit  value  is  found  which  reproduces  the 
experimentally  determined  But  this!  kind  of  search  increasofi 

imme-nsely  in  complexity  as  the  structure  of  the  medium  increases.  In 
order  to  use  this  method  pra.c  t real  ly , it  would  he  necessary  to  know  the 
*'i’  types  of  particles  in  the  mixture  to  he  modelled, 

as  well  as  the  for  all  the  types  hut  the  one  of  interest,  which  is 
the  value  that  this  method  would  he  nttem|)tin;’  to  delermiiu’.  Alterna- 
tively, and  much  more  simply,  ciiu’  could  i>erhai>s  expej  iment  al  ly  determine 

the  o^,'^  separately  for  all  the  censt  i tuents  of  a mixture,  and  then  the 
for  the  mixture.  From  (Al),  then,  for  an  ii-eomi)onent  mixture. 


f 

I 

I 


N. 

1 


dr  f.(r)  o^'^(r)  No 

1 ahs  ahs 


i/i  ■' 


I ahs 


(»8) 


This  eiiuation  has  a k.noun  value  of  n.  = N.o^,  ^ onlv  if  i)  it  was  directly 
‘ 1 .1  ahs 

mearuired  by  determin  Inp,  a.  in  a s i np,l  e-conponent  medium  (us.ually  not 
practicable),  or  il)  the  values  of  u . and  j / i,  were  di’terninrd 


experimentally  by  determininp  a and  , and  the  same  size  d 1st r ihnt ions 


( 

I 


id 


4 3 


were  kept  for  nil  e.oinponent  s t ]i  roup, limit . Then,  if  f^(r)  nncl  n.  nre 
kno  ,'n  ri‘ns;onnlily  well,  .in<l  k.  is  known  v.'ithin  or<i<'r-of-irap,ni  t iide,  tlie 
Mle  theory  p.oriLlnii.s  conlcl  he  operated  in  the  ,'ihovc-  .seareli  mode  to 
del  1‘rr.une  Tor  exampli-,  suppose  the  diffuse  refleelance  ('xperirienl.al 

test  uses  a sampT'’  whieli  is  a mixture  of  barium  suTpliate  and  lust  one 
other  component,  say,  carbon.  Imagine  that  o (carbon)  could  not  Ije 

«l  I)  s 

found  directly  by  experiment.  Then  the  first  equality  iu  Kq . (1)8)  would 

bo  used.  Nov,  If  component  number  1 is  carbon,  number  2 is  barium 

(2) 

sulphate,  and  if  = N„  o , and  ct  = No  , for  the  mixture  have  been 
2 2 abs  abs 

determined  by  reflectance  experiments,  tbcri,  from  Kq . (1)8) 


„(1)  ^ f 

abs  J ' 


dr  rpr)  - (»p 


■’[No  , 

abs  2 abs 


(1)9) 


and  the  search  mode  could  he  used  with  this  relation,  to  determini'  i 
for  carbon,  provided  f(r),  n,  and  N for  the  carbon  in  this  mixture 
were  quite  well  known.  If  all  the  carbon  particl  'S  had  (27ir/l)<<  1, 
then  the  first  term  on  the  TJIS  of  Eq.  (D5)  could  he  used,  v^^ei)',htod  with 


^1  ^ ■ 


where  v 


(1)  _ 

\abs r-  1)"<-27:)  . 

m +/ 


(1)10) 


dr  f (r)  4ti  r ^/3.  Th  i s 


are  kiiov.'ii. 


can  he  sol  veil  easily  for  k , if  n and 
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RADJATIVK  TKANSFKU  TiiKOKY  l.N  ONK  lUMKNSlON 

This  section  is  inclii(ii‘il  boc.iur.e  it  illustrates  simply  and  clearly 
three  very  important  fundamental  points:  i)  the  connection  between 

imaginary  refractive  index  and  absorption  cross-section,  ii)  the 
transition  from  v/ave  theory  to  radiative  transfer  theory,  iii)  the 
correct  difference  rather  than  differential  ecp-at ions  of  radiative  trans- 
fer theory. 

Let  the  strei  died  string  be  the  one-dimensional  physical  system 
which  supports  wave  motion.  Considin;  first  a uniform  string  of  mass/unit 
length  u,  under  tension  T.  This  string  supiiorts  transverse  oscillations 
u(x,t)  which  obey  the,  wave  equation 


(Kl) 


3t  9x 

and  vdiich  have  fundamental  solutions  n(x,t)  at  .a  given  free|uency  w: 


u(x,t)  = e^‘^*^(A^e  A_e^'^^) , k Y w/c,  c^  i/p,  (K?) 


where  A^  are  undetermined  compl ex- valued  amplitudes. 

Consider  anotlier  string,  also  under  tension  r,  but  witli  mass/unit 
leng.th  Pj  , with  oscillations  whidi  satisiy 


3^u  , / 3u  3^u 

"I  -T2*  '"d  ' ri  • 

3t  3x 


(F.3) 


Here,  the  term  /p^x  y 3u/3t  is  a frictional  term,  representing  absorption. 
The  fundamental  solutions  u(x,t)  at  f rispiency  m are 


• JT-f 
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I -imkx  . , imkx.  , , 

ii(x,t)  = e (A^e  + A_o  ),  k = oi/c,  (1 

with  m H /n-iyii/k  r n - Ir,  tc  > 0,  n = These  ehoicer.  nf  sij'.ii 

ensure  thnt  a traveling  wave  Is  clainpeil  in  its  direction  of  travel. 

A w;. ve  travel  in};  in  the  positive  x-direction  has  the  form 


iwt  -inkx  -2ttkx/X 
p e e 


(F.5) 


which  shows  that  the  wavespeed  on  this  string  is  c/n,  and  the  extinction 
coefficient  is  Attk/A.  For  small  Y.  (Vh/k  <<  1),  n ; n,  and  Amx/A  ; y* 

1.  Single  Scattering 


Now  imagine  a piece  of  the  frictional  string  of  length  r,  insertcal 
between  infinitely  long  segments  of  the  original  frictionless  string,. 

This  situation  is  exactly  an.ilog.ous  to  that  of  a particle  of  radius 

....  . , . / X iwt  -ikx 

r,  refractive  index  m,  in  vacuum.  het  a wave  Uj(x,t;  = A^  e e 

he  Incident  from  the  left  on  the  seg.ment,  .and  .i  wave  A^  e^*'* 

incident  from  the  right.  Boundary  conditions  ;ire  continuity  of  u and 

3u/9x  at  each  inter  f;ica-.  Application  of  these  boundary  condition.s  yields 

the  matrix  relations 


V2\ 


lA 


•f 


= (M) 


/ 


where  A^  is  the  coiaplex  amplitude  of  the  r i ght -t  ravel  ing. 


rig.ht  of  the-  segment,  and  A^  is  the  complc-x  .amplitude  of 


(I  f. ) 


w.ave  just  to  the 
the  Ic'ft- 


traveling,  w.ave  just  to  th(.  left  of  the  segment.  The  matrix  M i ;; 
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(M) 


1 -Hll 

1-n  I 

-i  kmr 

r 

0 1 /mil 

m- 1 

l-ri 

Ihii 

1 0 

ikmrl 

e ' \ m- 1 

ie'1 

0:7) 


Noti'  th.it: 

dot(K)  = ].  (E8) 


Consuler  the;  case  A2  = 0 (no  wave  incident  from  the  rij'.ht).  'J1u-n  Kq . 
(F.6)  yields 


R = I I I Aj^  I ^ = l^'*2 1 I I ^*22  ^ ^ ~ reflection  coefficient 

T r.  |a„|^/|a  1^'  = = transmission  coefficient 

(E9) 

R is  tlio  fr,  ction  of  the  incident  power  which  is  reflected;  T,  the 
fraction  which  is  transmitted.  If  m is  real,  R+T  = 1;  since  m is  complex, 
tlie  quantity 

ATI-  (R+T)  (nin) 

is  tlie  fraction  of  the  incident  power  which  is  ahsoriied.  Tliis  alisorption 
fraction  A is  tlie  .analop,  of  the  absorption  cross-section  of  a sc.it  t eri  np, 
particle  in  3 dimensions.  For  tlu“  case  k<<] , 2nr/A<<l , hut  n aiiiitrary, 
tlie  ahfivi'  equations  yield 


A = (4iucr/A)  (n)  . (i:i  1) 

For  n=l  , thi.s  is  the  exact  ;ina1o(’,  of  the  hy|>othesis,  (Fq.  D4),  ~ 

4tikv/A,  in  thrc'c  dimensions.  Here,  just  as  in  three  dimensions,  such  a 
hypothesis  is  badly  violated  in  p.eneral  . 
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2.  Trails  i L 1 on  to  Kncliative  Transfer  Theory 

Consider  a collection  of  the  ahsorhin;>  sefnix  iil  ;;  descritied  above, 
inserted  at  random  intervals,  vi.th  no  ovorlr'ip,  in  the  orij;ina1  fric- 
tionless string..  The  distr Unit  ion  c>f  these  sep, rents  is  to  he  honiO};en 
eons,  v.'hich  means  that  there  exists  an  averaj',e  i ntersep.ment  spacing,,  £ , 
and  thereby  an  average  number  of  segments  per  unit  length,  N 1/f, . 

Consider  the  n*^'*  segment.  bet  ^n+1^  amp!  i tr’es 

of  the  v/aves  incident  on  this  segment  from  the  (left,  rig.ht),  respect- 
ively, and  ^n+1^  complex  amplitudc-s  of  the  viaves  traveling 

away  from  (scattered  by)  the  segment  on  thc‘  (left,  right),  respectively. 

These  amplitudes  are  related  by  the  matrix  M,  as  in  Eti.  (Ff>).  Solution 
+ - + - 

of  Kq.  (K6)  for  (A  A ) in  terms  of  (A  , A , , ) yields 

n+1  n n nH 


/Cl' 

1 1 1 

1 1 

“21] 

^a"  1 

' n 

11 

b"2, 

, / 

(K12) 


wheri'  the  e(|uality  ^ = ~^'21  been  used.  These  amplitudes  A ..... 

are  the  amplitudes  of  the  various  waves  at  the  appropriate  interfaces 

betweiMi  tile  n*^'*  seg.ment  and  the  frici  ionless  string,. 

The  randomni'S!;  in  the  seg.ment  placing,  means  that  there  if;  zeri' 

correlat  ion  in  i)ha:;e  lietween  A , and  A . Tl'is  means  t'aal  •-  A A > --  (), 

' n 1 1 n ni  I n 


whefi'  < > expectation  over  all  posf;ilile  iiostions.  of  all  seg.ment! 

T.aking 


‘^|a  P>  1 . etc.,  and  using,  laps.  (K9)  , yields 

..  i.  1 » P f j * n ' n 


= Ti  *'  -i-  , 1 

nfl  II  n-H  n 


Ti~  . + Kl'’ 
nil  n 


(r.ria) 
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Solving  tlifsn  cqtiat.inns  for  tlu!  inLrnsit  i fs  I ' , yiold:: 

n"rJ. 


n+1 


n+] 


l'*'  = ^ t“  X -(l-T)l’*'  -t-  V.l  , 

11  n 1 n n n 


(Kl'J) 


<7  - f 


U-T)  I - KT', 
n n 


whom  the  approximate  equalities  follov;  for  R<<1,  A<<!  , T<1.  Here, 

the  intensities  T^.,,  I are  the  intensities  of  the  wave  disMir- 
n+1  n+1 

banco  at  any  point  botv%reen  the  n*'^'  and  (n+])*'^'  segment;  similnriy  for 

I*^,  I . This  follows  because  the  string  is  frictionless  between  seginents. 
n n 

Therefore,  if  the  mod  i un  is  iraaginc^d  to  be  divided  into  cells  of  length 

?,  = H \ with  on  the  average  one  segme.nt  anywhere  in  eacli  cell,  the 
+ 

iiiteiuiities  T ,,  can  be  taken  as  the  intinsities  on  the  right  cel]  houn- 
n+1 

+ 

dary,  v;I)ile  1^  are  tliose  on  the  loft  coll  boundary.  Dividing  both  sides 
of  Eqs.  (E13)  by  i,  treating  the  case  K<<1  , A<<1  only,  and  making  tlie 


transition  nP,  ->•  x - continuous  variatvle,  T~  ->  I j (x)  , yields 


dl 


■f 


dx 


(KiS)l^  + 8l_  , 


dl 

cfx 


(KlS)l 


ST 


4 ’ 


(in  4) 


wlieri- 


K KA,  S NIT, 


(in')) 


are  kept  constant  during,  the  limiting  process.  These  are  radiative 
transfer  tlieory  equations  for  Lliis  model  one-d  iinens.  i cnal  system.  Note 
that  they  have  the  same  form  as  the  KM  etiua t i our.  (Dl),  and  that  the 
s lg,n  i f icance  of  (K,S)  is  the  same;  In  liolh,  K --  aI);!or  pt  i<>n/un  i t lengtli. 


A 9 


S “ back-sc;if  tor  Lni /nil  i t length,  in  the  iroJium,  (cxo.opl  for  a fartor  of 
tvK)  which  is  geometrical  in  oiigln). 

Ilowovi.-r,  for  strongly  absorbing,  or  reflecting,  f:eg,nients,  tiie  full 
equalities  of  Kqs.  (hi  3)  nuist  be  used.  Letting  nf,  > x “ continuous 

-t  4 _ 

variable,  1 = I.(x),  1 ,,  - I^(x!-P),  these  eriiia:  ions  reduce  to 

n ± n+1  ± 

I^(x+i)  = (l-a)T^(x)  + bl_(x). 


I (x-ft)  = -bl^(x)  + (l+c)I  (x), 

— "T*  ~ 


(bib) 


wlierc 

a S (T-'f^+R^)/T,  b R/T,  c = (1-T)/T. 


(F.17) 


Equations  (Elf>)  are  difference  equations  rather  than  differential 

equations.  Tlie  desired  physical  solutions  are  tlie  "smoothest"  functions 

+ 

which  yield  the  correct  1^^  on  the  points  x ~ ni. 

These  equations  may  be  solved  by  normal  mode  metliod;.  li  t 

T)X 

I^(x)  = e'  in  E(|S.  (E14).  This  re.sults  in  the  eig.envaluo  equations 


(FIE) 


Solution  of  this  for  the  CMgenvalues  p,.  anti  c' i g,envectors  yield;: 


-1 


wlu're 


»'l,2  ^ ^-(1,;>)  ■■  ''  ^■'’^^i],7^\(l,3) 

1 1 ? 

<lj  2 ^ -^(f-a)  ' Ic-a  1 ^(c-a)^  . 


(FI  9) 
(FRO) 


The  fundamental  solutions  are  the  superpe;;  i t i on  of  the;.i'  normal  modes 
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p X p X 

1^(X)  - A_^,  <■  ' + c 

-1  I’]’' 

l_(x)  ^ b [(a+qj)A^j  n + (n-K|^)A^^  e ] 


(V.?]) 


It  is  slr.plo  nncl  str.i  i p,ht  for\v'rircl  to  vorify  that  these-  solutions  satisfy 
(the  first  c-qiinli  ties)  in  Kqs.  (K1  3) , on  the  iioints  x = n?. 

It  is  wc)rtli  noting  that  tlie  general  solution  of  I->|s.  (Klb)  may  he 
wr  i 1 1 cn 


_ e-i''  P^x 

I^(x)  = c fj^(x)  + e f2(x), 

-]  Pi''  P?'' 

I_(x)  = h [(n+qj^)e  "P  (a-lq2)e  ' f2(x)]. 


where  f^(x),  f2(x)  are  any  two  functions  periodic  in  x vjjth  period  C. . 

Slu'-.e  the  der.ircd  physical  solutions  are  to  be  the  "smoothest"  fvmrtions 

between  the  points  x = (nP.,  (n+3 )!’.),  it  is  necessary  to  choose  f^(x)  = 

A.,  = const.,  ItCx)  = A,„  = const.  The  possible  oscillations  represented 
"fl  / -rZ 

by  non-constant  f „(x)  result  from  tiie  extra  (infinitely  many)  di-grees 
J , z 

of  frc-c-dom  introduc(-d  in  going,  from  the  correct  discrete  ec|uations  (K1  3) 
to  tlie  continuous  difference  equation.-:  (I-ilb);  these  extra  degrees  of 
fri-edom  an-  clearly  spurious. 

It  is  instructive  to  note  that,  if  the  d i f f eii.-nce  equations  (l-ilfO 
wi-i'(-  fir.‘it  convertt-d  to  d i f fc-ix'nt  i a I equations  by  t iie  standard  method, 
us.ing,  hi  in(?.  >0)  1 ( i ^ (x+!  ) - T^Cx))/?]  dl^^/dx,  etc.,  and  then  sc'lved,  the 
results  would  iiavc-  be(-n 


dl  j^/dx  - f"'(-al^  I bl_),  dl_/dx  -=  y.-'(-hT^lcl_)  , 


(1-72) 


with  so1utir)ru; 
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q >•./?  q^x/P. 

I^(x)  =■■  e ^ ■ 


_i  qi^/^- 

T_(x)  - b l(a+q^)A^j  c + (a+ci2)A,^e  " ] 


(r,2  3) 


wliere  Liri(£  >0)  [n/P , h/t,  c/f,]  = constants.  Insfxct  ion  of  I'.qs.  (ri9), 

(K21)  sliows  tliat  tlic'so  expressions  are  tlie  correct  solutions  only  if 

hf  since  tiien  ^ ^ ?.n(Hqj  ~ q^  ,^/i..  Tliat  is  to  5;ay, 

only  if  the  clianp.cs  in  over  one  cell  arc'  small  are  the  differential 

equations  (E23)  equivalent  to  the  difference  equations  (T16).  In  order 

for  these  chanp,es  to  he  small,  |q  | must  he  <<1  . Tf  both  A<<]  , !!<<!, 

1 , / 

then  the  differential  equations  (I’23)  are  indeei!  just  tlie  equations 
(ElA).  Note  that  |q^  2^^^^  follows  if  only  A<<! , re};ardless  of  R.  if 
only  A<<]  , l)ut  R < 1,  the  differential  equatiems  (R22)  are  correct,  hut 
they  arc‘  not  tlie  same  <as  the  Ktl  e(|uations  (I'l4).  This  ))oint  is  made 
clear  by  an  example  in  Appendix  1. 

The  solutions  (1021)  of  tlie  difftrenee  equations  (Klh),  subject  to 
the  boundary  conditions  analap.ous  to  tho.se  for  a plane  i>ar;illel  s.imple 
of  Lhichness  d in  3-dimeiisions, 


l_(c1)  = 0,  1^(0)  = 1 


(H2^0 


yield  the  followinp.  expressions  involviu!;  the  ref  I ect  ancc's  R^l , R^_ 

1-1’  .i:,,, 

R^  = (al-q^)/!.;  (p,-!'^)d 

d “ 


O’ 2 3) 


\/luTo  K,  1 (0)  f(>r  a sair.i>lo  ol  total  d . In  tbo  limit  A<<!  , 

(1 

H<<!  , ther;o  reduce  exactly  to  tlie  form  of  the  KM  la  lations  (U^)  , 


r-  - 


I 
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I 


K 

S 


(1-n 

Of> 


K 


1-RjK 

d 


Sd  n 

1 - 1>  d 


(F26) 


CO 

where  K,S  nre  defined  liere  by  Kqs.  (r.l5),  wirli  N - <1“^. 

It  is  enliy.htenlnj;  to  contrast  the  tv;o  solutions  liy  means  of  a 
numerical  example.  Consider  strongly  reflecting,  moderately  absorbing 
segments,  with 


A = 0.1,  K = 0.8,  T = 0.1.  (E27) 

Then,  from  Fqs.  (i;l7),  (f'.19),  (E20), 

a = 7.3,  b = 8,  c = 9;  ^ 2.40,  -0.70; 

2 = 1-22/P-,  -1.20/?..  (K28) 


The  correct  relations,  from  Fq.  (F25) , are  then 


R 

Co 


1-R  R 

0.825,  ?.n(-- - 2.4?  d/?. 


(R29) 


11  the  standard  KM  relations  (F20)  were  used  with  tlu’se  numerical  v.tlu(':;, 
the  rel.'itions  (F25)  would  read 

1-R  R 

R^^  = 0.01,  = 0.8?'.  d/f  (KIO) 

which  are  s i g.n  i f 1 cant  1 y in  c>rror.  (>n  the  oihei  1>  u’d  , if  tljc  incorrect 

.solutions  (F23)  wt'ie  used  in  the  honnd.iry  v.iliu’  problein,  without  re<iuirin,, 

|q,  R would  come  out  to  he  the  cot  root  v.ilue  (lK’‘t),  since'  it 

1,2  “> 

floes  not  flepend  on  lu  . However,  the  second  rel.it  ion  in  (F.25)  would  be 

L , / 


(1 


1.1  d/f. 
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diffc'Vf’nt  from  both  tlie  iiocoiid  rolntion.s  in  (K29)  and  in  (F30)  . 

In  sum,  in  lb  i;.  appendix  it  lias  boon  shown  that,  lor  a ono-climcns  i on,-i] 
random  scattorinp  modium, 

1)  Kadiat  Lvo  transfer  (-(luations  follow  directly  from  wave  t. henry 
equations  if  there  is  zero  phase  correlation  amonp,  the  waves  incident 
on  any  scattcrer  from  different  directions. 

li)  Tlie  correct  radiative  transfer  equations  ;ire  difference  ratlier 
than  differential  equations. 

This  latter  point  is  exploited  in  Appendix  I. 
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MTKNniX  F 

CONVKRSJ.ON  OF  FRACTION  BY  VJFICHT  TO  NUMBFK  OHNSIIY 

In  a nulLi-confJt  i tuoiiL  random  scattc-rinf.  medium,  tlie  (;iveri  data 
are  invariably  tbe  total  mass  density,  I),  of  i lie  actual  medium,  and 
the  fraction  by  wcifjit,  of  each  constituent  number  i.  For  use  in 

Mie  theory  alp.orithms,  it  is  necessary  to  V;now  the  number  density 
of  each  of  the  constituent  ])articles. 

Let  be  the  mass  density  of  particles  of  type  j,  and  be  the 
average  volume  of  a type  i particle.  Then,  for  a sample  with  total 
mass  M,  total  volume  V,  mass  of  i*^'^  constitui-nt , and  n^  “ total 
number  of  type  particlc’S, 


m 


1 

V 


II  .P  .V. 
3 1 1 

V 


N .p  . V . . 
Ill 


(FI) 


Solvlnp,  this  for  yields 


N.  = ii.D/p.v.. 

1 3 11 


(F2) 


Now 


V . 
1 


dr  f ^ (r)  (4TTr^/3) , 


(FI) 


where  f.(r) 

3 

in  order  to 


is  the  s i x.c"  distriluit  ion  of  type  i jiart  icle.';. 
evaluate  , the  size  distributions  must  also 


Therc'fore , 
be  iven. 


5>5 


APPI'NDTX  G 
SIZK  niSTR]  Bin' IONS 


Four  woll-known  particle  size  dist  rihut  ion<-  were  proj’r.'inimed  for 
use  In  tills  work.  These  .are  the  two-par.iineter  pamiaa  d i st  rlhution , the 
bltnodal  gamma  distribution,  the  log  normal  distribution,  and  ,a  power 
lav;  distribution.  ITie  unnomualized  distributions,  n(r),  and  the  nor- 
malized ones,  f(r),  are  defined  and  described  below. 

1)  Two-parameter  gamma  distribution 


n(r)  = 


I 1/b  -3  . / , V 

|r  r exp-(r/ab),  r^^  < r < r2 


0 


, r < r^  ; r>r^ 


If  r.  = 


- 0,  r^  = '”,  then 


(Gl) 


a = 3<’r>  - 2r  , b - o^/(l+2o^),  - (<r^>  - <r>^)/<r>^ 

o 

where  r - mod.al  radius  = a - 3ab , <r>  = mean  radius  = <a-2.ab,  <r^>  = 
o 

2 1 1 

mean  square  radius  = (ah)  (j^  “ ^ ^ ~ *"^>0,  or  both, 

then  the  relations  between  (<r>,  and  (a,b)  are  altered,  but  r = 

a-3ab  unler.s  r,  > r or  r„  < r . 

1 o 2 o 


ii)  liimod.il  ',.Tnn’:.i  distribution 


n(r)  = 


n(r ,a^ ,b) 


n(r,a,,  ,b) 


rr 


— 4 


(G2) 


dr  n(r,aj^,b) 


dr  n ( r , ;i 2 , b ) 


1 1 

where  n(r,<n.,b)  is  :t  two-par;iPiet  er  gamma  distribution  with  jiarameters 


(.a^,b).  Here,  h.'ilf  of  the  particles  are  in  eacli  mode. 
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111)  Lo{’  iiorni  il  cli.sti  ihnt  ion 


-1 


r (!xpri-(^n  vh'^)  r^<ri'r^ 


n(r) 


(C3) 


r<r^,  r>r^, 


If  = 0,  = ",  then 


^ ^ -0^/2  p 2 -2o^ 

<r>  = r e , <r^>  ~ r v , r 


r e =-  <r 

{; 


O'  . p^l/2 


and  |(<r^>  - <r>^)/<r>^|  ^ J-e  ” . If  r > 0,  or  r„  < ",  or  both. 


these  rrl.ifions  are  altered. 


iv)  Power  lav;  distribution 


fr  , r . < r < r 

‘ run  ~ max 


n(r)  =«i 

lO  , r < r . , r > r 

mn  max 

2-a  2-a 

1 t-  - r . 

, ^ ^ I -a  max  min  j , o 

here  <r>  ~ — — , a / 1,2. 

2-a  1-a  l-a  • 

r - r 
max  min 


(CA) 


<'r>  = 


2-a 

r 

max 


2-a 
r . 
ni  I n 


(2-a)en(r  /r  . ) 
max  mm 


- , a = 1 


(l-a)S’.n(r  „ /r  . ) 


<r>  = 


max  min 


l-a  l-a 
r - r . 
max  min 


, a = 2. 


The  normalized  distributions  f(r)  wi're  ealeulateil  numeiic.al  ly  from 

To 


f(r)  - n(r)/ 


dr  n(r) . 


(fl5) 


.Several  other  distrihut  ions  could  be  .addtd  easily  to  tlie  rom]'uter  proj;ram. 


Arrr.f<mx  i> 


bi 


I'UOM  p(ro.-.O) 

111  tlris  v/ork,  tlic  sinp.lo  sr.it  tet  i ii};  ph.isi'  fuiulioii  p(cf>sO)  v,;is 
fir£.'t  c.'ilciil.ntcd  from  !Uc  tlu'ory  for  so.  :o  set  of  sr.il  tiT  iii>',  .iiiplo:;  , 
k = 1,.  kTi.it  is  iierdc'cl  for  ir.iiltiplr  sr;it  tor  i np.  r.il  rul.ni  ions  in 

til  is  study  is  the  tiz  iniutli-.-ivcr.'ipoil  iih.nse  function 

2ir 

) ~ ~2~  I ‘^'1''  p(cosO),  (HI) 

o 

on  some  sc-l  (p^.  i,.i  = 1,...,  N.  Hero, 


cosO  = Pit'  + /j-p^  *'^-p'^  cosiji'.  (H2) 

The  notliod  used  in  this  work  to  ohtfiin  p^^^^(p,ti')  is  tlmt  of  Ikinson^  ^ 

iind  Iknnsen  .'ind  Tr.svis*.  First,  choose  the  p.ir  t i ciil  .nr  (p.,p'.)  v.nlues 

1 :i 

desired.  For  e.ich  p.nir  of  v.nlues,  do  the  i ntepr.nt  ion  in  I'ci.  (Ill) 
nun.er  ic.nl  ly.  For  cx.ninple,  iisinp,  C.nuss-Lep.rnd I'e  (|u.idr.nt  ure  on  t he 
points  (^'  “ <1^,  with  weights  v;  , n"l,...,Q, 


(p.,t>'.)  tv  p(p.p'.  + /l-tl?  COSlIi  ) 


(113) 


However,  the  v.tlu(>  cor.'':  ..  = p ii ' -f  >^l-ii'.^  co.sA  nip.ht  not  be 

n 1 .1  i I i j n 

one  of  the  rosO,  for  wliicli  p(cosO  ) was  c.n  1 cii  1 .it  eil . Tf  not,  1 inee.r 
K K 

interpol.nl  ion  is  used  hetwetm  the  two  .ndiacent  cosO,  to  find  i>(cos'^  ..). 
' k ' nil 

After  p^"^(Pj,p'.)  is  found  by  this  method  on  .nil  points  (vij,ti'.) 


which  .nre  desired,  it  is  r enorm.n  1 i zed  liy  usinp, 
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.1 

)/|  8p  p^*’\p  ,P  ' ) . (il‘'0 

-1 

The  function  p^°\)i,p')  so  defined  is  thus  correctly  nornu'il  i x.(td  to  w^. 
This  ronormnlizati on  helps  to  remove  er»-ors  flint  may  have  hoen  intro- 
duced by  the  interpolation. 

This  method  of  obtaininj’  p^'*^(fi,p*)  is  considerably  faster 
computationally  than  the  method  wliich  uses  Legendre  function  expansions, 
for  a given  accuracy. 


I 


.'j9 

ait);m)tx  j 

NKi;  .K(il!AT10N  OF  KADTAXm,  TKANSFI'R 

In  th<-  iisiinl  (Icr  i v.'it  ion^  of  tlu'  sLnnd.ird  RT  ftiiintion  (A9) , tlicrn 
is  no  nnaloj;  to  tlic  cUfforonce  ec]uation.s  (F16).  Howevf'r,  th'-  derivation 
of  F{q.  (A9)  given  in  Appendix  A docs  pass  througli  suci)  an  analog,  namely, 
Eq.  (A6).  The  considerations  in  Api)endi>;  E lead  to  tlie  conclusion  that 
the  diflcrcitce  equation  (A6)  ought  to  he  a more  accurate  radiative 
transfer  equation  than  the  standard  differential  equation  (AO),  for 
close-packed  strongly  absorbing  or  strongly  reflecting  media.  Therefore, 
(Ah)  is  tentatively  considered  as  the  new  equation  of  radiative  trans- 
fer. 

In  what  follows,  consider  only  prohltmis  which  arc  invariant  to 
translations  nornal  to  the  x-direction.  Then,  the  intensities  in  Eq. 

(A6)  may  be  written  1 (zl?.p/?,p,4’)  , etc.,  whcie  p is  the  cosine  of  tlte. 
angle  between  s and  the  x-di  rect  ion,  and  <})  is  tlie  azimutlial  angle  of  r. 
around  the  z-;ixis.  For  applications  in  vdiich  only  the  tot.al  tran;;mitted 
and  reflected  fluxes  are  needed,  the  (})-depcndenco  of  Eqs.  Af>  may  be 
integrated  out,  just  as  v/as  done  in  going  f ror:  I’q.  (Alh)  to  Eq.  (A18). 

Eq.  (A6)  t lien  reduces  to 

1 

I(z-Kp/:i,p)  - (l-ye)!  (z-£p/2,n)KYi/2)  dp’  p^'’\p,p’)l(z-Cp’/2,p’). 

-1 

(11) 

Since  z is  a cont  inuous  variatile,  let  z > z + lf.li/?  in  Et].  (II). 


The  resulting  equal  ion  is 
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] 


I(z+£ii,u)  = (l-Y£)I(^-,,ii)  + (y£/2) 


dp'  P^"^P,p’)r((zf  yP(p-p'),!-'). 

(T?) 


■1 


Let  there  be  nn  incident  beam  in  direction  y.  at  z - 0,  f;o  that 

o 


I.  (o,p)  = ^ F fi(p-p  ).  (13) 

in  2 o 

Tiien,  let  satisfy  I(2+P,y,i0  = (1-yP)  T(z,i0.  Tlie.n 

I.  (z.p)  = -•  F i5(p-p  ) exp[z  r(u)],  (lA) 

in  2 ® 

where 

I'(n)  = (p£)  £>.i(]-y£);  (13) 

note  that,  for  » F(p)  ~ ~ y/p>  the  usual  form.  Then,  let 


l(z,p)  = T.  (z,p)  + I (z,p)  (rf.) 

in  d 1 

in  Kcj.  (12),  where  T^.(z,)i)  is  the  diffuse  inttnsity.  'lliis  yields 

Iji(x.+£u,p)  = (1-y"  ) j ('^  »p)  + (y*-  F/A)>/°\i!  ,\i^^)  exp[r  (v^)  (?+  -|l.(p-Vi^))  1 

1 

+ (y£/2)|  dll'  p^°^  (p  ,p' )T^I^  (z  + if  (p-li' ) ,li' ) . (17) 

-1 

This  differi'uce  ecjuat  ion  is  the  analnj;  of  the  d i f f e riMi t i a 1 eiiuation  (Al'f). 

It  is  not  clear  tu'W  to  solvi'  this  equal  i('u  hy  a dcudiliiu;  metliod, 
wliitli  depends  on  tlie  ehoiee  of  a start  Inp,  layi-r  so  thin  tlial  tiie  last 
term  ahovi>  may  he  droppi'd  for  the-  transfer  across  that  layer  (scc‘ 

Appi'iidix  J).  Here,  the  miniiiRiiii  thickness  layer  which  can  he  chosen  in 
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principlf  is  one  of  tliickncss  9.,  bt'caiisc  that  is  tlio  tliickness  of  a 
moiiopnrt  i cl  c-  liiycr,  and  tl:  v/holc  derivation  of  tlic  oquntions  in  terms 

of  parLiclc  cross-sections  probably  nabos  :•<>  sense  for  thinner  layers. 
Now,  if  which  implies  that  both  a <<  1 , B <<  1,  then  tlie  differ- 

ence etiuation  (17)  reduces  to  the  standard  differential  equation  (A?1). 
But,  if  accuracy  of  1%  or  better  is  desired,  then  the  conventional 
equation  (A21)  .sliould  not  be  used  for  y!.  > 0.01. 

A many-flux  solution,  similar  to  that  done  by  Mud';ett  and 
Klchards^^  for  the  standard  diffetential  eqiuition,  may  certainly  bo 
done;  this  is  quite  cumbersome  and  timt — consuming  for  more  than  four 
fluxes.  However,  a combination  of  an  accurate  many-flux  (16  or  '12-flux 
might  be  sufiicient)  solution  for  the?  starting  layer  of  thickness  Z, 

VJith  a doubling  method  solution  for  the  larger  thicknesses  desired, 
seervs  feasible. 

A Monte  Carlo  solution  is  of  course  possible.  In  fact,  it  seem?: 
likely  that  the-  standani  Monte  Carlo  technic]ue  is  already  solving,  this 
difference  equation,  rather  tlian  the  usual  differential  eejuation, 
inasmuch  as  the  Monte  Carlo  technique  looks  at  tlie  radiative  transfer 
collision  by  collision.  If  so,  discrepancies  should  exist  between 
Monte  Carlo  metliod  ifsults  and  results  of  other  metiiods  for  the  st.indaril 
equat  ion,  for  closely  packed,  riedi.a  cont  .ai  n i ng,  strong.ly  absorbing,  or 
strong,ly  re  f I c'c  t i ng.  p.ir  t i c 1 es. . This  pcjint  rem.iins  to  be  checked  in 
future  work. 

Tile  two-flux  t heoi  y of  K<|uation  (12)  may  be  obt.ained  in  a 
manner  similar  to  th.it  in  wliich  two-ilux  theory  of  liq . (A1B),  tlu‘  KM 
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tlu'ory,  v;ns  obtaiiicil  in  Appcnolix  I’>.  Hnko  tlio  ;ipprox  iri.il  ion.‘ 


T(z+P.p,|  ) - I,  (z+Pp)  - I_^(zff./2),  p > 0 


I(z-Pp.-p)  - ] (z-Pp)  = T (z-P/2),p  > 0 


vhero  thp,  necorKi  equnlit  ios  follow  from  tlio  dcfini.t  Lon  p jclii  V'  “ • 

o 

Then  Fq.  (12)  reduces  to  the  pair  of  coupled  equations 


I^(z+P/2)  = [I-(aFp_)P]l^(z)  + (F_P)  r_(zt-P/2)  , 


I_(z)  = (P_P)l^(z)  + |l-(a+r._)P]]_(z+P/2),  (TO) 


± ■ 2 


-1  1 

fdp  fdp'  p^"\p,'p’). 


just  as  in  Appendix  15,  and  a,  fi,  y have  the  same  rusininns  as  in  Appendix 
A,  with  (5,  + [5  = P,  so  = a IP  . 

T"  ~ 'r  ~ 

Define,  as  in  Appendix  15, 


.1. 

't’^(z)  E I dp  P I,^(z)  = Y (111) 


'Jtten  15(is.  (10)  reil\ice  to 


‘*>|(z+P/2)  = T 'lj(z)  + U <!_(z-L'P/2) 


^ K '1,  (z)  + T <T_(zl  f/2) 


(112) 


T " l-(ntp  )f,  K p P,  A . c«P, 


6'i 


satisfy  A + K f T = 1,  and  tiavf  tlic  siy,i\i  f ic.a\icc‘  (T,  1?,  A)  = (tran.:- 
niissieni,  reflection,  al'scirpt  ion)  fractions  for  a layer  of  thicknc.ss  P./2. 
Tlie  presence  of  H./2  r.at.lier  titan  H in  the  arpnnent  s of  is  a result 
of  the  averap/inp  over  ii  done  in  ohtaininj;  Kr|s.  (T9).  Note  tliat,  if 
R <<  1,  A <•'  1,  then  Kqs.  (T12)  reduce  to  the  KM  ccpiations  (Rl),  as 
cxpectf'd;  the  presence  of  X./2  ratlier  than  ?.  in  (112)  is  crucial  for  this. 

Solving  these  eqtiations  for  <>^(z+t/2)  in  terms  of  'i>^(z)  yields 

4'.(z+!l/2)  = (l-a)0,(z)+h  <t  (y.) 

+ + - 


(II  4) 


<!■  (z+i/2)  = -h  0,(z)  + (l+c)0  (y) 
— + — 


wlier  o 

a ^ (T-T^+R^)/T,  b r R/T,  c - (1-T)/T,  (11 ‘>) 


just  as  'n  Kqs.  (K17).  Since  thette  diffc’renci'  eciuations  (114) 
same  forri  ;is  the  Rqs.  (I'lh),  the  results  and  solutions  of  Rqs. 
may  be  applied  directly  here.  llte  solutions  are 

PjZ 

\(y-)  - e h 

r-/ 

<’_(/.)  = b [(a+(|p  A^l  e + (alq^)A^^  e |, 


h.ive  the 

(in  4) 


(rih) 


where  A,,,  A,_  arc-  tiwlet  ermined 
1 J r/ 


coc  f f i c i cut  s. , 


and  t.'hei  i 


>’1,2  ' '‘l,2^’ 

9l  2=1  ' K'-'-)  (117) 


for  a si. lb  of  tliickiK 


d ('ll  ;i  total  l\  absculi  i ng,  I'ackg.ronnd , 


M 


1-K  R 

Thf.se  are  just  the  rc^RuJ!.,  !:lh-i:2]  and  K25,  cxc  ;'t  for  tlit:  factor  of  2 
in  the  definitions  of  2‘  1 f T < ] , R <^  ] , A <<  i , i'>is.  (118)  ;»o 

ever  to  the  KM  Theory  expressions  (B2) . 

A fev;  cxar.iples  are  mo.st  instructive  here,  and  slun,'  tint,  in  some 
pra.ctica]  lu'ohlems,  the  results  predicted  hy  the  corrected  equation.s 
(nA,  118)  differ  significantly  from  those  predicted  by  the  stamlard 
KM  equations  (B1 , R2) . 

i)  Consider  a powder  made  up  of  small  carbon  particles.  For 
such  particles,  R <<  1 , but  A < 1,  so  T<<  1.  Let  T ' 1-A-R,  so,  from 
Kq.  (115), 


c-a 


2 2 
(A+Il^  t-  Jl 

1-A-R 


_jA  _ RA(2-A) 

(1-A)’ 


Ol'^) 


to  0(R).  Then,  for  A ~ 1/2,  say,  c-a  is  not  <<1  , so  q is  c|uit(' 

^ y ^ 

diffeii-nt  from  (c-a)  , and  p / 2q  This  ii;eans  that  the 

I j / J-  > ^ 

rel.'ition';  (T18),  among  coef  f i c i t'lit  s (8_,a),  are  q ite 

different  from  the  KM  relations  (R2).  A value  of  A > 1/2  is  quite 
realistic  for  c.arbon  powdi-r. 

ii)  ('.onsider  the  riodil  jiuvders  t ts.‘a.t  i d in  t lii  s work.  In  all  of 
thorn',  t'  ■ ;.tron;;l'.'  ri' f 1 ec  t i ng.  I'.afO^  war;  bv  far  the  na'M'rity  conq)f'nent  . 
For  example,  in  :;amplc  11,  lla.fit'^  t carlion,  the  valuer;  of  (p  ,n,C)  may 
he  found  f ri'in  Tabler;  2,  1,  at  wavelength  1 = 0.6  pm: 


n - ) . / 


1x10  (imi) 


0.  102 (pm) 
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There foro,  from  F,q.<;.  (113), 


K 0.283,  A = .513x10  T = l-R-A. 

Since  A <<  1,  we  expand  the  parameters  (a,b,c)  of  Kqs.  (115)  to  0(A) 
The  results  are 


R . A(1-2R+2K  ) 


(1-R) 


2 


K 

(1-R) 


„ ^ .2RA.1/2 

,,2  ’ 'll, 2 ■■  ■ ^1-R^  l’l,2  ""  ^■'ll,2^^' 


Tlien,  from  Kqs.  (112) 


')u(^  R'>  1/2  1/2  1-R.R 

_ .2AU-R),.  4d,2RA,  „ , d » , 

R ) .-TW  = 

a CO 


Solving  these  for  A/R  and  R yields 


1-R  ^1  4cl  R 

('ir)A  - —2 — 1 (it 


(120) 

Co  (1  oo 


Tliese  relations  are  to  be  contrasted  with  Eqs.  (B2)  of  the  KM  theory, 
which,  for  A<<1,  yield 


A ^ K 
R S 


(1-R  )‘ 


lid  1 1 

-T  ^/l^,^/R  )■ 

CO  (j  CO 


(121) 


Note  that  A = a9.  = lCe/2,  R = p_J  = St/2.  Comparison  of  (120,  121) 
makes  It  clear  that,  for  given  data  (R,,R  ), 

Cl  <o 


(t — g'l  ~ (ll)i/M  ; or  •'  l’>0 

1-R  new  KM  new  RM 


<")new  = (^>KM- 


I 
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That  i*',  thr  KM  thc-ory  prediction  of  K in  too  lrij',li,  Iiy  a factor 
l/(l-0.28'i)  ~ 1.4  in  this  numerical  or.ampio,  Init  ttio  KM  tlioory  prodiction 
of  A ap.rees  witli  tlint  of  the  correctod  tiioory.  'J'hi;;  latter  oqvinlity 
i.s  fortunate  for  the  main  purpose  of  this  v.'orl;,  wliich  vat;  to  check  tiie 
supposed  equality  of  tlie  KM  coefficient  K = 2A/f.,  calculated  from  the 
KM  eqs.  (It?),  wi  tli  the  absorption  coefficient  2a.  Note  well  tliat  even 
thoup.li  |q  9 iti  this  cxampl  e,  tlie  results  do  not  agree  with  those  of 

the  KM  theory;  for  this  agreement  it  must  he  true  that  both  A<<1,  R<<1, 
which  was  not  the  case  in  cither  of  these  examplet!. 

It  is  worth  remarking  that,  for  very  strongly  reflecting  particle's, 
especially  for  such  particles  of  la  ge  radius,  the  usual  assumption  of 
neglecting  the  mnltiiile  scattering  in  a nionoparticl.e  layer  may  be  in- 
valid. 'Ihis  assumption  was  implicit  in  the  above  derivation  of  the 
corrected  two-flux  equations.  In  order  to  take  into  account  such  effects, 
at  le.'ist  ;i  1-flux  set  of  equations  must  be  used,  or,  for  accuracy,  a. 
many-flux  set.  It  is  also  likely  that  the  conventional  reasoning  which, 
allows  replacement  of  all  the  different  particles  in  the  medium  by 
average  particles  (see  Appendix  A)  needs  to  be  reexamined  if  y!  •>■  1. 

These  are  topics  for  future  work. 

In  .';um,  the  above  two-flux  theory  of  t lie  difference  equation'. 

(12),  and  the  two-flux  theory  (>f  the  d i f f c'l'ent  i a 1 eipiation  (A18),  (the 
KM  theory),  yield  significantly  dilfi-rc'nt  results,  for  strongly  absorb- 
ing or  .stron;;ly  reflecting,  closily  packed  ii;edia.  This  implies  that  the 
accurate  solutions,  of  ia|.  (17)  and  liq.  (Al'l)  will  be  s i g.n  i f i cant  ly 


different,  for  finch  media. 
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Prcvimts  n t ! f-i;ip  1 5;  l>y  otlir>r  .nui  hors  at  discontinuous  thoorios  of 
radiative  I -ansfiT  are  outlined  in  reference  (11).  These  llic-ories  are 
only  the  two-flux  approximations,  and  the  individual  particle  .scatterinp, 
and  absorption  cross-sections  do  not  seem  to  be  an  intep.ral  part  of 
the  development,  although  too  philosojiby  of  these  theories  is  similar 
to  that  of  the  new  theory  prc.sented  here. 


I 

I 
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Al’f'KNmX  J 
nOUBT.INf;  HKTIIOn 


In  this  apjKMul  ix,  the  douhlinj’  incaliod  c-quation.s  used  in  Llils  work 
arc  drrivt‘(!.  Tlie  f inidanental  transmission  equations  for  a plane  para- 
llel layer  of  tliicknesr.  f,  (optical  thickness  t = yH)  may  be  written 

1 

= [dll' fT^(n,u')I^(o,n')  + Rj^(ii,n')I_(d,ii')] 
o 

(.11) 

1 

l_(o,u)  = |dp'[Tj(li,,i’)I_(d,p')  + k^(p,ii')I^(o,i.')l 
o 

where  )j  > 0,  and 

I^(z,p)  r 4i  (z,p) . (.1?) 


Here,  I^(z,p)  arc  the  total  intensities,  diffuse  plus  directly  t r.ins- 
mitted,  averap.ed  over  azimuthal  anple  (J , and  are  defined  hy 


T^(z,l0  = I(z,p),  p>n;  I_(z,ii)  ? l(z,-u),  ii>n.  (.11) 

In  symbolic  operator  form,  F(|s.  (.11)  imiy  he  written 

I,a)  = 1^(0)  + I_(D 

(.IM 

(_(o)  = Tj,  I_(0  i Rj,  T^(o). 


Consider  another  layer  of  thickness  d,  between  z = f and  z - ?£.  Then 


= T|,  I^(e)  t-  R^  7_(:>1>) 
7_(P)  = Tjj  t R^  7^(£). 


(.l''>) 
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But,  for  the  total  tlilc’cncss  29, 

= '^29.  '^29 

(J6) 

L<o)  = ^29.  ^-(29.)  + I^(o). 

Combininj^  these  sets  of  equations  yields 


which  are  the  doubling,  equations  in  operator  form.  I’or  example. 


X 

- Id,l' 


1 

dp"  T^(P’,P") (...), 


(J7) 


v.’horc  the  ellijisis  (...)  stands  for  the.  opt;rand,  some  function  of  ii" . 
The  operator  ^ defined  by  its  series  expansion. 


(1-R  R ) 

9.  l' 


-1 


ItR  R„  d R R R R\  + ... 

9,  9 9 9 9 9 


(.18) 


jd,.*  fi(p-p’)(...)  + dp'  Rj,(p,p')|dp"  Ry(p',p")(...)  + 


In  pr.u'fical  ce'rrput  at  i on , ten  or  fewer  terns  ari'  kept  in  tiu'  above- 
series,  and  an  estimate  of  round-off  c^rror  is  includeil,  followinp.  the 
method  of  liansen  and  Travis*. 

In  oilier  to  make  use  of  the  doublinj’  equations,  the  t ransf(>r 
m-itriees  , R^  for  ;i  very  thin  layer  of  thieknes;;  9,  so  thin  that 
multiple  scatterinj’,  m.iy  be  nej'.lec  t ed , mur:t  be  found.  I'or  such  a l.iyer, 
the  scatteriiif;  term,  “7|dp'  p ^^'^  (|) , p ' ) 1^1  . (z  , p ' ) , is  dropjied  from 
I’.qs.  (A19),  and  the  resulting,  equation  is  solved  analytically.  llie 
T9,  Ry  may  be  read  off  from  the.sc'  solutions;  they  turn  out  to  be 


I 
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i I 


Tlu-si*  c:-:pr('.S5;i  c>ns  iirc*  finite  for  all  cons  idi’ialil  c‘  care  miisf 

l>e  taken  in  computation  ns  eititer  p or  p'  ->■  0,  or  as  )i  >■  p'. 

Numerical  i ntc-p,ra  t f on  over  (S-functlon  i nti'j'.raiul!!  can  lend  to 
sij'.nif leant  error.  It  wan  judp.od  more  accurate  to  write 


^£,d  i ^ 


(.110) 


at  each  step  of  the  douhllnf',  where  contains  the  fi-fuiiction  oper- 

ator- This  eximndn  the  eciiiatlon  for  to  fovir  teins,  and  that  for 
to  five.  For  doulilinj’  from  a thickness  I.  to  a thickness  ?h,  t'nero 


resul ts 


hu  - T,^  - r, . 


(dll) 


where  T„,  , . , T,  contain  all  the  terms  not  proportional  to  f(p-p'). 

21,,  di  l,,di 

If  thi'  Incident  intensity  is  <^(p-p  ),  then  th('  incident  flux  in 

o 

p . For  a naniple  of  thickness  d,  the  diffuse  transmittial  and  refli'Cted 
o 


fluxes  arc 


I 

i.,.  <l|i4i  (,(D,p)  ^ dp/ii  T ,,()i,p  )>'Ti 

+(!  i r J (1  , (!  t O (> 

o 

I 

o_d-  = jdpv'i;  i_(i).p)  - jdp/;;  '^ct,<ii^“’"o^’^o 


(M2) 


7] 
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AJTKNDfX  K 
COMPUiLK  PROCitAMS 

Two  conpu'.er  prof.rnnis  were  developed  duriii;'  the  course  of  this 
investlf’ation , one  based  on  the  Hie  theory,  tlie  othc'r  on  conventional 
multiple  scatterlnf;  theory.  A short  description  of  each  program  is  given 
in  this  appendix.  Complete  documentation  and  listings  are  available 
from  the  author. 

Both  programs  use  Gauss-Legendre  quadrature  and  double  precision 
aritlimetic  throughout,  and  both  were  thoroughly  checked  for  accuracy. 

Tlie  Mie  scattering  program  was  checked  against  published  results,  and 
against  analytic  results,  and  it  was  verified  that,  ns  the  model  sample 
particles  are  taken  further  and  further  apart,  the  near  field  phase 
function  (Appendix  A)  approaches  the  far  field  one.  The  multiple  scat- 
tering program  agrees  with  the  results  of  an  analytic  calculation^^,  and 
satisfies  internal  cliecks  for  all  optical  dtq'ths  < 1024. 

1 . Mie  Scattering  Program 

Given  tlie  mass  density  of  a model  mixture  of  spiiericnl  particles, 
and  the  optical  constants,  size  distriljution,  and  fraction  by  wclglit 
of  each  component  in  tlie  mixture,  tills  program  calculates  the  nuraiuT 
density  of  tlie  mixture;  tlie  mimber  density,  near  and  far  field  differential 
cross-section,  and  the  extinction,  scattering,  and  alisorption  cross- 
sections  for  each  component  of  the  mixture,  and  for  tlie  mixture  as  a wliole, 
using  standard  Mie  tlieory  a1  gori  t inns , and  tlie  j)rocc>dures  explained  in 
Appendices  A and  F.  It  also  calculates  the  piiase  functions,  on  a cliosen  set 
of  scattering  angles  0. 


73 


A dctailrd  dcscri]>tion  of  a typical  Mic  tlumry  scattorinp,  c.'.lcula- 

3 A 

tlon  Is  readily  available  ia  the  current  literature  ’ . Tlu;  major 

di.f fcTcnces  between  the  proj;ram  developed  for  this  work  and  the  Pr.AUSf) 

17 

routine  are  as  follows: 

i)  The  logarithmic  derivatives  of  i|/^(z),  rather  than  the 

functions  themselves,  wore  computed  by  hackxi^ard  recursion  ^ ^ ^ Here, 

the  are  the  Riccati-F.essel  functions,  the  notation  that  of  Van 

3 

de  Ilulst  • 

11)  Tlie  output  of  the  program  yields  the  phase  fi:nction  p(0)  on  a 
selected  sot  of  scattering  angles  0,  not  the  coefficients  of  the  Legendre 
polynomials  in  the  standard  series  expansion  of  the  phase  function. 

2.  Mviltiple  Scattering  Program 

This  prog.ram  uses  the  punched  output  (mixture  i)hase  function, 
p(0))  from  the  Mic  scattering  program,  and  computes  the  diffuse  reflec- 
tance and  transmittance  as  functions  of  the  angle  0 of  the  incident 

o 

beam  and  the  plane-parallel  sample  optical  thickness,  using,  the  doubling 
method  (Appendix  J)  to  solve  the  standard  radiative  transfer  equation 
(A21).  The  program  then  uses  the.'-e  calculated  reflected  fluxes  in  the 
KM  formulas  (II?)  in  order  to  compare  directly  the  predicted  values 
(K,S)  with  the  defined  values  K r 2h,  S “ 2(1_.  The  program  also 
calculates  the  quantities  (P^)  > defined  in  Appendices  15  and  C. 

The  phase  function  p^”^(p,p')  is  obtained  from  the  input  p(0)  by  tin' 
method  described  in  Appendix  H. 


r 


